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Abstract 

We provide an algebraic classification of all supersymmetric domain wall solutions of 
maximal gauged supergravity in four and seven dimensions, in the presence of non- 



trivial scalar fields in the coset SL(8, JR)/SO(8) and SL(5, JR)/SO(5) respectively. These 
solutions satisfy first-order equations, which can be obtained using the method of Bo- 
gomol'nyi. From an eleven-dimensional point of view they correspond to various con- 
tinuous distributions of M2- and M5-branes. The Christoffel-Schwarz transformation 
and the uniformization of the associated algebraic curves are used in order to determine 
the Schrodinger potential for the scalar and graviton fluctuations on the corresponding 
backgrounds. In many cases we explicitly solve the Schrodinger problem by employing 
techniques of supersymmetric quantum mechanics. The analysis is parallel to the con- 
struction of domain walls of five-dimensional gauged supergravity, with scalar fields in the 
coset SL(6, lR)/SO(6), using algebraic curves or continuous distributions of D3-branes 
in ten dimensions. In seven dimensions, in particular, our classification of domain walls 
is complete for the full scalar sector of gauged supergravity. We also discuss some general 
aspects of D-dimensional gravity coupled to scalar fields in the coset SL(N, H)/ SO(N). 



1 Introduction 



Recent years have seen an increasing interest in gauged and ungauged supergravities 
in various dimensions following the conjectured duality between gauge theories and 
string/M-theory |TJ, |^, ^J. The AdS/CFT correspondence offers the possibility to un- 
derstand strongly coupled gauge theories from a dual supergravity description. The 
relevant backgrounds have Poincare invariance along the boundary directions and are 
asymptotic to Anti-de Sitter (AdS) space. The deviations from the AdS geometry in 
the interior correspond either to the broken phase of the theory due to non-zero vacuum 
expectation values of scalar fields 0], [§J-fl2| or to deformations of the conformal field 



theory 13|-[21 



In this paper we shall be concerned only with the first possibility. In particular, we 
will construct a large class of solutions of D = 7 and D = 4 supergravity that are dual to 
the (2, 0) theories in six dimensions |22| and the three-dimensional theories with sixteen 
supersymmetries |23| on the Coulomb branch. These theories correspond to the world 
volume theories of parallel M5- and M2-branes respectively. In contrast to the conformal 
cases, the solutions include also non-zero scalar fields, which, although they vanish at the 
boundary, become large in the interior. These solutions can be lifted to eleven dimensions 
describing the gravitational field of distributions of a large number of M2- or M5-branes. 
The location of the branes is directly related to scalar Higgs expectation values on the 
field theory side. Using such solutions, the spectrum of scalar and graviton excitations as 
well as the expectation values of Wilson loops can be calculated, shedding new light on 
the AdS/CFT correspondence. Such investigations have been carried out for supergravity 
duals of Af = 4 SYM, sometimes with surprising findings || |9|, [TIJ. Interestingly, many 
of these backgrounds arise as limits of charged AdS black holes in the lower-dimensional 
gauged supergravity theories or as extremal limits f)| |6|, PH of rotating brane solutions 
24] , [|, |25, ; they have null singularities near the continuous distributions in the higher 



dimensional backgrounds. Such singularities arise generically in the flows from confor- 
mal to non-conformal theories and it would be interesting to understand them better 
from the field theory side.Q Furthermore, these geometries can be viewed as examples 
of consistent truncations in various dimensions. The embedding of lower-dimensional 
into higher- dimensional supergravities has been worked out explicitly only for a couple of 
cases including the S* 4 and S 7 compactifications of eleven-dimensional supergravity p8| 
to D = 7 |29|] and D = 4 |3(J gauged supergravity, respectively. More examples of explicit 
Kaluza-Klein Ansatze which relate gauged supergravities to ten/eleven-dimensional su- 



pergravity have been worked out in the recent papers [31 



1 Recently, supergravity duals of supersymmetric gauge theories with eight supersymmetries were 
constructed and a mechanism for resolving a certain type of space-time singularity was proposed [p7| . 
However, the situation in 27 is different because, in TV = 2 SYM, the classical superconformal point at 
the origin of the Coulomb branch is removed by quantum effects, whereas the Coulomb branch of M = 4 
SYM is uncorrected. Therefore, it is not expected that a similar mechanism is at work for backgrounds 
with sixteen supersymmetries. 
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An additional motivation for investigating solutions of gauged supergravity in detail, 
as well as the spectra associated to quantum fluctuations, stems from the possibility 
to apply them to scenarios that view our world as a membrane embedded non-trivially 
in a higher dimensional non-factorizable space-time |32| . This old idea has been revived 



recently in relation to the mass hierarchy problem (see also ||34|| ). In the latter work a 
slice of the AdS$ space, where our four dimensional world is embedded, was cut out, thus 
resulting in a normalizable graviton zero mode. However, it turns out that there exists 
in addition a continuum of massive graviton modes with no mass gap separating them 
from the massless one. It was shown in [^] that these modes have negligible effect to 
Newton's law. However, it is quite desirable in this context to find ways to model possible 
modifications of Newton's law, since any deviations from it at the sub-millimeter scale 
have not been ruled out by the present day experiments (see, for instance, f35j). It is 
known, on the other hand, that there are Yukawa-type modifications to Newton's law 
in theories that possess a mass gap separating massless from massive graviton modes. 
Besides the obvious phenomenological advantage of such models, the existence of a mass 
gap leads to a well-defined effective field theory of the standard model of particles plus the 
massless graviton. This mass gap should be independent of the details of the slicing or at 
least be practically insensitive to them. Toy models possessing such desired features have 
already been constructed in |I6] in the context of five-dimensional gauged supergravity. 
In this sense, many of the models, that we will describe in the present paper, as well as 
many of the models in [[□]], can also be used to further pursue these ideas. 

This paper is organized as follows: In section 2 we describe the bosonic sector of 
gauged supergravity in D dimensions. The non-zero scalars take values in the coset 
space SL(N,]R) / SO(N) common to gauged supergravities in any dimension. We re- 
strict to vacua with (D — l)-dimensional Poincare invariance and find that for certain 
values of D and N the equations can be cast in first-order form. These Bogomol'nyi- 
type equations are equivalent to Killing spinor equations and these solutions preserve 
sixteen supersymmetries. In section 3 we present solutions in arbitrary dimensions that 
preserve part of the SO(N) isometries and examine some of their general features. We 
explain how these solutions, for the cases that correspond to gauged supergravity, can 
be lifted to string/M-theory, thus showing that they are consistent truncations of ten- 
or eleven-dimensional supergravity respectively. These higher- dimensional backgrounds 
arise in the AdS/CFT correspondence as supergravity duals of the field theories living 
on D3, M2 or M5 branes on the Coulomb branch. In section 4 we give a general dis- 
cussion of the scalar field equations and graviton fluctuations in these backgrounds, and 
the mass spectrum of operators in the dual field theory. Our discussion is facilitated by 
the fact that the corresponding equations can be cast into one-dimensional Schrodinger 
equations with appropriately chosen potentials. We also make contact with the theory of 
supersymmetric quantum mechanics and outline the necessary elements that will be used 
for solving explicitly the Schrodinger problem in some cases of current interest. Section 
5 is devoted to a detailed discussion of the differential equation for the conformal factor 
of the background metrics written in conformally flat form. We present a complete clas- 
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sification of solutions with different unbroken isometries in terms of irreducible algebraic 
curves. The problem of solving the domain-wall equations is then essentially reduced to 
the uniformization of the relevant curves and the inversion of the corrsponding functions. 
In Section 6 we treat in detail distributions of M2-branes, which correspond to curves 
with genus g < 1. We present their uniformization and obtain explicit expressions for 
the conformal factor and the Schrodinger potential of the equation for the scalar and 
graviton field fluctuations, whenever this is possible in analytic form. In almost all cases 
of genus zero, the spectrum can be found exactly, whereas for the other models we uti- 
lize the WKB approximation to determine the mass gap and spacing in the spectra. In 
section 7 we repeat the same analysis for distributions of M5-branes. We also include, as 
application, the calculation of vacuum expectation values of Wilson surface operators in 
the six-dimensional (0, 2) theories on the Coulomb branch using the eleven- dimensional 
backgrounds. Section 8 contains a summary of the Lame equation and its generaliza- 
tions, which arise in the study of scalar and graviton fluctuations in the background 
of domain walls associated to elliptic functions. Elements of supersymmetric quantum 
mechanics are also used to expose connections between different elliptic potentials, and 
discuss briefly some features of their exact spectrum beyond the WKB approximation. 
Finally, we end with a short discussion in section 9 and list some open problems. 

The present paper generalizes previous work by two of the authors fLl| , where all 
domain-wall solutions of five- dimensional gauged supergravity with non-trivial scalar 
fields in the coset SL(6, 1R) / SO (6) were classified in terms of algebraic curves. In various 
places we include for completeness, but with no further explanation, results on the D = 
5, N = 6 theory, which correspond to continuous distributions of D3-branes in ten- 
dimensional type-IIB supergravity. 



2 Gauged supergravity and first-order equations 

We consider the sector of gauged supergravity theories with non-trivial fields in the coset 
space SL(N, 1R)/ SO(N). This is in general only a subset of a larger coset space and is 
common to gauged supergravity theories in any dimension. To be specific we will study 
the following Lagrangian in D dimensions 



1 1 N ^ 

C = - A n--Y,{da I f-P{a I ) , (2.i; 



i=i 



although we will be mainly interested in the two cases (D, N) = (4, 8) and (D, N) = (7, 5) 
that correspond to a sector of the four- and seven-dimensional gauged supergravities 
|30| , p9[ . Another case of interest, corresponding to a sector of the five-dimensional 
gauged supergravity |37|, j38H , is (D,N) = (5,6). All other fields including the fermions 



are zero. In this subsector the scalar potential P can be expressed in terms of a symmetric 
N x N matrix M = S T S, where S is an element of SL(N, H). Using an SO(N) rotation 
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matrix, M can be diagonalized 

N 



M = diag{e 2/ \...,e 2/3jv } , ^# = 0. (2.2) 



In addition, since the determinant of M is 1, it depends on iV — 1 independent fields. We 
parametrize the fields $ by the N — 1 independent scalars aj in the following way 



N-l 



(3i = Y^ x ^ a i , (2-3) 



i=i 

where Xn is an N x (N — 1) matrix. The rows of this matrix correspond to iV weights of 
the fundamental representation of SL(N), which obey the following normalizations: 

N-i 2 N N 

i=i i=i i=i 

Then the potential takes the form 

P = -£ [(trM) 2 - 2tr(M 2 )] . (2.5) 

The equations of motion that follow from varying the action (|2.1|) with respect to the 
metric Gmn and the scalars ai are 

1 1 N ^ 1 

-TZmn — - dM&idNOti — — — -GmnP = , 

i=i 

ftp 

d M (^GG MN d Nai ) - V^Gj— = . (2.6) 

ooli 

For the applications we have in mind we need the metric to exhibit (D — l)-dimensional 
Poincare invariance, namely 

ds 2 = e 2A(z) {ri^dx»dx u + dz 2 ) = 

= e 2A{r) ri^dx»dx u + dr 2 , (2.7) 

where the ralation between the coordinates z and r is such that dr = —e A dz. Furthermore 
we assume that the scalars depend only on z and we concentrate on solutions arising from 
first order equations. 

The traditional way would be to consider the Killing spinor first-order equations in 
the three different cases, when ( |2.1D corresponds to a sector of gauged supergravity, as 
it was done for (D,N) = (5,6) in ||. That would leave half of the supersymmetries 
unbroken. However, there is an alternative way to arrive at first-order field equations 
using a method a la Bogomol'nyi f39|-||42|l that will be useful for further generalizations. 
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For this task we have to plug the ansatz (|2.7|) into the action and rewrite the potential 
in terms of the prepotential 



4 

We find 



1 1 N 

\V = --trM= -jJ2 e2Pi ■ ( 2 - 



The resulting one-dimensional effective action 

S = / dre ^(i:(^) 2 ->-l)( D -2)g) 2 + 2 p), (2,0) 
can be written as a functional, which is a complete square plus a boundary term 




+ 2g(D-2)e< D - 1 > A W\~_ x , (2.11) 

provided that the number of scalar fields and the dimension of space-time are related as 

D — 2 

N = 4- — -. (2.12) 
D-3 K J 

This relation will be used to simplify various expressions in the rest of the paper. 
Then, from ( [2.11| ) we can read off the first-order differential equations: 

dA q doir q dW 

^ = -D^2 W ' ~dr ~ \~dai ' ^ 



Note that ( |2.12| ) has integer solutions only for the values of N and D that were mentioned 
at the beginning of this section, namely (D,N) = (5,6), (D,N) = (7,5) and (D,N) = 
(4,8). In these three cases there exists a maximally supersymmetric solution of the 
equations of motion ( |2.6|) that preserves all 32 supercharges, in which all scalar fields are 
set to zero and the metric is just AdSo- In these cases the potential in (|2~5| ) becomes 
P = — g 2 N(N — 2)/32 and equals by definition to the negative cosmological constant. 
This defines the mass scale g that we have already used. The associated length scale R 
defined by g = 2/R will also be used in this paper. The same AdSr, space is obtained for 
general values of D and N when all scalars are set to zero, but then, there is no notion 
of supersymmetry. 



3 The general solution 

We begin this section with the construction of the most general solution of the non- 
linear system of equations within the D- dimensional ansatz ( |2.7|) that preserves Poincare 
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invariance in the embedded (D — l)-dimensional space-time. To keep our considerations as 
general as possible we will assume that the scalar field potential is given by (|3.1|), but we 
will not assume from the very beginning that the relation between the number of scalars 
and the dimensionality of the space-time (|2.12|) holds. After discussing some general 
properties of the corresponding configurations, we concentrate on those cases where the 
relation (|2.12|) holds, and therefore the solutions originate from gauged supergravity in 
the appropriate number of dimensions. Then, we lift our solutions to eleven or ten 
dimensions in the context of eleven- or ten-dimensional supergravity. Our analysis is 
based on an analogous treatment of the case (D, N) = (5, 6) that was developed in 
reference fill. 



3.1 Solutions in arbitrary dimensions 



It is possible to find the most general solution of the coupled system of equations ( |2.13| ). 



As we have already mentioned, we will consider the potential arising from expression 
1]D after using ( pT2|) . We find that 



Unless otherwise specified, we will no longer restrict N and D to obey ( 2.12|) . The 
resulting theory is that of iV — 1 scalars coupled to D-dimensional gravity with the 
interaction potential ( |3.1|) . However, even though it no longer represents, for arbitrary 
values of N and D, a sector of some gauged supergravity theory, it still admits interesting 
solutions arising from the first-order Bogomol'nyi equations ( |2.13| ). 



In order to proceed further, we first compute the evolution of the auxiliary scalar 
fields pi. Using fl2j) and (^TJ) we find 



ft = 2^A> + 9 -e^ +A , i = 1, 2, . . . , N , (3.2) 

where the prime denotes derivative with respect to the argument z. It is easy to inte- 
grate these N decoupled first-order equations for the /Vs. For further convenience we 
reparametrize the function A(z) in terms of an auxiliary function F(zg 2 ) as follows 

e A = g(-F') hd-v+n ? (3.3) 

where the prime denotes here derivative with respect to the argument zg 2 . The minus 
sign we have included in this definition implies that, for consistency, the function F 
should be monotonously decreasing with z. Then, according to this ansatz we find that 
the general solution for the D-dimensional metric ( |2.7| ) is 

ds 2 = g 2 f^^r] IMU dx"dx u + g- 2 f-^dF 2 , (3.4) 

with 

N 

i=i 
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and the solution for the scalar fields in ( |3.2| ) is given by 

fl/N 

e 2ft = ^-— , i = l,2,...,JV. (3.6) 

r — Oi 

The bi's are N constants of integration, which can be ordered as 

h > b 2 > . . . > b N , (3.7) 

without loss of generality. Also, since the sum of the is zero, we find that the function 
F has to satisfy the differential equation 

(-F') A = TT(F - h) = f , A = < D -V N . (3.8) 
11 ' A(D-2) + N 



If we assume that our models arise from gauged supergravity, and therefore the rela- 
tion (|2.12|) between the number of scalar fields N and the dimensionality of space-time 
D holds, the metric (|3.4|) will become 



ds 2 = g 2 f^^r] liu dx fl dx u + g^f^^dF 2 , (3.9) 



whereas the scalar fields are still given by ( |3.6|) . The function F(zg ) is determined by 
solving the equation 



iv 



(^ = 11(^-6,) = /, (3.10) 

i=i 

since the constant A equals 4 in all three cases (D, N) = (4, 8), (5, 6) and (7, 5). Although 
the relation (|2.12j) guarantees that our models originate from gauged supergravity, it 
should be emphasized that it is not the only possibility that results into an integer 
value for A. As an example, choosing D = 5 and N = 4 we obtain A = 3. Another 
interesting case, but with non-integer A, arises for D = 3 and N = 2: it has A = 4/3 
and describes three-dimensional gravity coupled to a single scalar in the presense of a 
negative cosmological constant (since the potential (|3.1| ) turns out to be constant). 

So far, we have presented our general solution in a coordinate system where F is 
viewed as the independent variable.0 If we insist on presenting the solution in a confor- 
mally flat form, as given by the first line in ( |2.7|) , the differential equation ( |3.8|) needs to 
be solved to obtain F(zg 2 ). This will be studied in detail in section 5, using the theory 
of algebraic curves and their uniformization, as it is a necessary step for investigating the 
Schrodinger equations that arise for the massless scalar and graviton fluctuations. 

If the constants 6j are all equal, our solution becomes nothing but AdSo (with radius 
2(D — 2)R/N) with all scalar fields turned off to zero. In the opposite case, when all 
constants 6, are different from one another, there is no continuous subgroup of SO(N) 
preserved by our solution. If we let some of the 6,'s coincide, we restore various continuous 

2 To compare this with the results of |fl| for the values D — 5 and N — 6, we should replace the 
function F(z/R 2 ) and the constants 6, used in [O by AF(zg 2 ) and Abi, respectively. 
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subgroups of SO(N) accordingly. Imposing the reality condition on the scalars in (J3T 
restricts the values of F to be larger that the maximum of the constants bf, then, according 
to the ordering in ( |3.7|) , this means that F > b\. For F ^> b% the scalars tend to 
zero and / ~ F N , in which case the metric in (|2.7| ) approaches AdSn (with radius 
2(D — 2)R/N) as expected. In the conformally flat form of the metric, where z is viewed 
as the independent variable, we have e A ~ 1/z and therefore z = is taken as the origin 
of the z-coordinate. Hence, in solving the differential equation ( |3.10| ), we will choose the 



constant of integration so that in the limit F ^> b\ we have the asymptotic behaviour 
F ~ l/z 4 ^^ 2 )/^.^ For intermediate values of F we have a flow in the D-dimensional 
space spanned by all scalar fields In general we may have b\ — b 2 — ■ ■ ■ — b n , with 
n < N, when b\ is n-fold degenerate. In this case, the solution preserves an SO(n) 
subgroup of SO(N) and the flow is actually taking place in N — n dimensions. On the 
other hand, let us consider the case when F approaches its lower value b\. Then, the 
scalars in (|3.6|) are approaching 



fi /N (F-hY n ~ N y N , for i = l,2, 



n 



e 2ft ~ { } , (3.11) 

,1/N 

{^(F - b x ) n ' N , for i = n + l,...,N 



where f = n^n+i(^i — Consequently, we have a one-dimensional flow in this limit 
since the scalar fields $ can be expressed in terms of a single (canonically normalized) 
scalar a as 



a , y^plHF- bl) . (3.12) 



When our solutions correspond to gauged supergravities and the relation (|2.12|) ap- 
plies, it is also useful to find the limiting form of the metric (|3.9| ) when F — > b\. Changing 
the variable to p as 

N 

F = b 1+ ([l~n/N)fl /N gp) N - n , (3.13) 
the metric (|3.9|) becomes, for p — > + 



ds 2 ~ dp 2 + ((1 - n/N) n f g {D - 3){N - n)+n ) t^- 3 ^-) pW=m^r) ViV dx»dx v . (3.14) 

Hence, at p = (or equivalently at F = &i) there is a naked singularity that can be 
interpreted, as we will see later in a higher-dimensional context, as the location of a 
distribution of M5-, M2- or D3-branes for (D,N) = (7,5), (4,8) and (5,6) respectively. 
It will also be seen that this singularity is time-like for n — 1, 2, 3 and null for 4 < n < N. 

3 Using we can show that the scalar fields go to zero as ai ~ z D ~ 3 , when ( 2.12 ) is obeyed. This 



is consistent with the interpretation that these scalars parametrize states in the Coulomb branch of the 
N = 4 SYM theory and they do not correspond to explicit mass deformations. 



S 



3.2 M-theory branes 



It is possible to lift our seven- and four- dimensional solutions with metric and scalars 
given by ( |3.9j ) and ( |3.6| ) with (D,N) = (7,5) and (4,8), respectively, to supersymmetric 
solutions of eleven-dimensional supergravity, where only the metric and the three-form 
are turned on. The eleven-dimensional solutions will correspond to the gravitational 
field of a large number of M5-branes and M2-branes in the field theory limit with a 
special continuous distribution of branes in the transverse to the brane space. It shows 
that they are true compactifications of eleven-dimensional supergravity on S A or S 7 , 



respectively; this is also what is expected on general grounds |43| . The method we follow 
is the same as that used in [Tl[] to lift the five-dimensional solution with (D, N) = (5, 6) 
to a supersymmetric solution of the ten-dimensional type-IIB supergravity representing 
distributions of D3-branes in the field theory limit. A brief summary of the main results 
will also be included for completeness in this case. 

The higher-dimensional metrics for the various distributions of branes have the form 



M5 — brane : 



ds 2 



H lf \ v dx^dx v + H^{dyi + dy 2 + . . . + dyg) 



r 2/3 



M2-brane : ds 2 = Ho 2/ \ u dx^dx u + H^(dy 2 + dy 2 



1/3/ 



dy\ 



(3.15) 
(3.16) 



and 



D3 — brane 



ds 2 = H l,2 r]^dx^dx u + E^\dy{ + dy%+... + dy 2 ) . (3.17) 



r l/2 



In all cases Hq is a harmonic function in the iV-dimensional space JR N transverse to 
the brane parametrized by the yi coordinates. However, instead of being asymptotically 
flat, the metrics ( |3.15|) - (|3.17|) will become asymptotically AdSp X S N ~ X for large radial 
distances, with D and N taking their appropriate values. The radius of the sphere 
is always R, whereas, in agreement with our previous normalization, it is \{D — 3)R 
for AdSr>. Under these conditions, the higher-dimensional solutions break half of the 
maximum of 32 supersymmetries (see, for instance, [Q). It is only for coinciding branes, 
when the metric is exactly AdSp x S 1 ^ -1 , that the maximal number of supersymmetries 
is preserved and the backgrounds are presumably exact vacua. 

We proceed further by first performing the coordinate change 



Vi = 2g 



(D-5)/2 



{F-bt 



,1/2 



1,2,. ..,N , 



(3.18) 



where the x^s define a unit iV-sphere. It can be shown that the iV-dimensional flat metric 
in the transverse part of the brane metric ( |3.15| )- (|3.171 ) can be written as 



N 



2 _ D -6 



N „ 2 N 

Y,^dF 2 + Ag D -^F-h)dx 2 . 



i=i 



i=l 



(3.19) 



i=l 



The harmonic function Hq is determined by comparing the massless scalar equation for 
the eleven- and ten-dimensional metrics ( |3.15| )-( |3.17| ), with the same equation arising 
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using the D-dimensional metric ( |3.9|) . In both cases one makes the ansatz that the 
solution does not depend on the sphere coordinates, i.e. $ = e lk ' x <j)(z). Since the solutions 
for the scalar $ should be the same in any consistent truncation of the theory, the 
resulting second-order ordinary differential equations should be identical. A comparison 
allows to determine the function Hq as follows 



N 



The coordinate F is determined in terms of the transverse coordinates yi as a solution 
of the algebraic equation 

N 2 

£^- = v 5 - ( 3 - 21 ) 

i=l 

Note that the algebraic equation ( |3.21| ) for F cannot be solved analytically for general 
choices of the constants \. However, this becomes possible when some of the b^s coincide 
in such a way that the degree of the corresponding algebraic equations is reduced to 4 or 
less. We also note that Hq as defined in ( |3.20| ) and ( |3.21| ) is indeed a harmonic function 
m U N . The proof was made in [|TT| for the case N = 6, but even if N is kept general all 
intermediate steps are essentially the same. We note that one may use as independent 
variables, instead of the y^s, spherical coordinates since the constraint ( |3.21| ) is then 
automatically satisfied. For D3-branes, this was done for the various cases of interest in 
|PT| . For M2- and M5-branes we expect to recover the metrics found in ||12|| . Various 
subcases were also considered before in || H, |8|, (| in connection with the Coulomb branch 
of gauge theories at strong coupling. 

Brane solutions that are asymptotically flat are obtained by replacing H in (|3.15 )- 
( p,17| ) by H — 1 + H . Then, in this context, the length parameter R has a microscopic 
interpretation using the eleven- dimensional Planck scale lp or the string scale ^/a* and 
the string coupling g s , and the (large) number of branes iVb. For M5-branes we have 
R 3 = 7rAy|, for M2-branes R 6 = 32nN h l P , and for D3-branes R A = 47rN h g s a' 2 . 



4 The spectrum of fluctuations 

In this section we investigate the problem of solving the differential equations that arise 
for the massless scalar field as well as for the graviton fluctuations in our general D- 
dimensional background metrics (|3.9| ). After presenting some general features, we make 
contact with supersymmetric quantum mechanics, which will prove useful for making 
further progress in computing the exact spectrum later in sections 6, 7 and 8. 



4.1 Generalities 



Within the AdS/CFT correspondence |], ||, §] (assuming ( 2.12|) ), the solutions and eigen- 
values of the massless scalar equation have been associated, on the gauge theory side, 
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with the spectrum of the operator TrF 2 , whereas those of the graviton fluctuations po- 
larized in the directions parallel to the brane, with the energy momentum tensor T^ v 
|4"5| , 0, |3|, P> . A priori these two spectra are different. However, it can be shown, sim- 



ilarly to the case of five-dimensional gauged supergravity investigated in [p6|, that the 
two spectra and the corresponding eigenfunctions coincide. Hence, in what follows, $ 
will denote either a massless scalar field or any component of the graviton tensor field. 

We proceed further by making the following ansatz for the solution 

z) = exp(ik ■ x) exp ' 

which represents plane waves propagating along the (D — 2)-brane with an amplitude 
function that is ^-dependent. The mass-square is defined as M 2 = —k ■ k. Then the 
equation for ^ can be cast into a time-independent Schrodinger equation for a wave 
function ^(z) as 

- ^" + = M 2 ^ , (4.2) 

with a potential given by 

v JD-tf 2 D-2 

4 2 v ; 

It should be noted and will be further discussed later in this section, that this potential 
is of the form that appears in supersymmetric quantum mechanics and therefore the 
spectrum is non-negative [47|, f5|.[] However, this does not guarantee that the massless 



mode is normalizable. In our case, \l/o ~ exp((£) — 2)A/2) and it is clearly not square- 
normalizable due to its behaviour near z = 0, namely \l/o ~ 1/ z^ " 2 ^ 2 . Using the general 
formulae (|3.3|) and (|3.8|) it turns out that the potential takes the form 



f2/A 

V = 

4i? 4 



N * „ N 



' i=i i=i v 



,2 



(4.4) 



We note that in deriving (|4.4|) from ( [4.3|) the relation ( |2. 12|) has not been used. 

The expression ( |4.4j ) for the potential depends, of course, on the variable z through 
the function F(zg 2 ). Even without any knowledge of the explicit ^-dependence of the 
potential, we may deduce some general properties about the spectrum in the various cases 
of interest. In general, F takes values between the maximum of the constants bi (which 
according to the ordering made in ( |3.7|) is taken to be hi) and +oo. When F — > +oo (or 
equivalently z — > + ), the space approaches AdSo and the potential becomes 



D(D-2) 1 

V ~ — i— , as z -> 0+ , 4.5 



and hence it is unbounded from above. Let us now consider the behaviour of the potential 
close to the other end, namely when F — *■ b\. Consider the general case where b\ appears 



4 In the context of five-dimensional gauged supergravity this has been first hinted in ]Ii[ | and explicitly 
noted in lEl. 



11 



n times, as in the corresponding discussion made at the end of subsection 3.1. Using 
( |4.4| ) we find that the potential behaves (including the subscript n to distinguish among 
various cases) as 



,2/A 

V n ~ J -<±— + A)n 2 -%n){F- b^ 2 , as F-6i, (4.6) 



with /o being a constant given, as before, by fo = Yli(bi — h). Hence, for A < n < N, 
the potential goes to zero and the spectrum is continuous. For n = A the potential 
approaches a constant value, which is given by VA, m in = 4]=^ A 2 /q^ A . Therefore, although 
the spectrum is continuous, it does not start from zero, but there is a mass gap whose 
squared value is given by the minimum of the potential. Hence, for A < n < N the range 
of z necessarily extends to +00, i.e. < z < 00, with z = 00 corresponding to a null 
naked singularity (except in the AdS case where n = N). If A < n < N we have that 
(F — &i) 1-n / A ~ (n/A — l)g 2 f^ A z near F — b\, and therefore the potential behaves as 

C 

A < n < N : V n ~ , as z — > 00 , 

f l/n 2 

For n = A we have instead that F — b\ ~ e - -* 9 z , as z — > 00 with the potential reaching 
the constant value that we have mentioned above. For n < A the potential goes to 
either +00 or —00, as F — ► b\ and therefore the spectrum cannot be continuous but 
discrete. Therefore there should be a maximum value for z, denoted by z max , that is 
determined by solving the algebraic equation F(z ma _ x g 2 ) = b x . The value of z — z max 
corresponds to a time-like naked singularity. Then, using the relation (F — bi) 1 ^ 11 ^ ~ 
(1 — n/ A)g 2 f^ A (z mSiX — z) near F = b\, we find that 

n< A: V n ~ — , as z -> 2 max , 

'max J 







- ^ 


1 


/ n 


64 





" 64^ A^ J "4- (48) 

One may worry that there are cases where the potential goes to —00 at z = z max with 
a coefficient C n ^ that is smaller than —1/4; then, it is well known from elementary 
quantum mechanics, that the spectrum might be unbounded from below. However, 
this does not happen because we are dealing with supersymmetric quantum mechanics 
and the spectrum is bounded from below by zero. In addition, the coefficient C n ^ in 
( |4.8|) is manifestly greater or equal than —1/4; the limiting value —1/4 is reached for 
n = 4 A/ (A + 4). The qualitative analysis of the spectrum we just presented agrees with 
that performed for (D, N) = (5, 6) in and for (D, N) = (7, 5) and (D, N) = (4, 8) in 

As we see later, solving the Schrodinger equation (|4.2j) and determining the spectrum 



is a non-trivial problem, except for a few particular cases where explicit calculations can 
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be carried out in detail. However, when the spectrum is discrete, as for n < A, we 
may use the approximate method of WKB. Since the quantum mechanical potentials are 
super symmetric, WKB is expected to be an excellent approximation, not only for high 
quantum numbers, but also for low ones 13]. Moreover, if the potential turns out to be 



shape-invariant, the WKB approximation is also exact ||48|| . It is convenient at this point 
to use F as an independent variable instead of z. Then, changing the dependent variable 
in (|4.2j) as \l/ = e^ D_2 ^ A / 2 0, we find the equation 

g 4 d F f l ^d F (j) + M 2 f^-^(f) = . (4.9) 

Using well-developed methods for studying this type of differential equations (see, for 
instance, fl25|), the spectrum is found to be given approximately by 



9 7T 2 ( D-3 l(n/4- l)A + ra| . „ , 
M 2 m = —m[m+^-+ 11 7 2(A _ n) ') 0(m°). = 1.2 (1.10) 



where, as usual, z max is the maximum value of z and A is the constant defined in ( |3.8| ). 
It turns out that the validity of the WKB approximation requires that the inequality 
n < A be satisfied. This is also consistent with the fact that the potential should be 
unbounded at the end points F = +oo and F — b\. 



4.2 Relation to supersymmetric quantum mechanics 

We have already mentioned that the potential (|4.3|) has the form encountered in super- 



symmetric quantum mechanics [47, 48|. Let us make this relationship more explicit by 
first recalling that in supersymmetric quantum mechanics, two potentials are supersym- 
metric partners of one another provided that there is a superpotential W(z) so that 

Vi(z) = W 2 - W , V 2 (z) = W 2 + W. (4.11) 

Then, in terms of the ladder operators 

a = ^- + W(z), a f = --!L + W{z) , (4.12) 
az az 

the eigenstates of the two Hamiltonians Hi, H 2 are related to each other as 

i-S> = -r^r^Si . €1 = ~^="Hi? , (4.13) 

whereas for the energy levels, in the case that the spectra of the two partner potentials 
are discrete, we have the relation 

EW=E£li, n = 0,1,2,..., (4.14) 

with Eq = 0. Notice that if V'i+i °^ Hi ^ s normalized then the wave-function tpn^ of H 2 
will be also normalized and vice-versa. Thus, on general grounds, knowing the spectrum 
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of Hi one can construct the spectrum of H 2 ; conversely, from the spectrum of H 2 one can 
deduce the spectrum of Hi apart from the ground state with = which is not paired. 
This relation is true only for the case of unbroken supersymmetry. When supersymmetry 
is broken, there is a 1-1 pairing of all eigenstates of Hi and H 2 and the relations become 



modified by replacing tp^-li whh ipn } and E£j.i with En' in the equations above. Then, 
the potentials V\ and V 2 have degenerate positive ground state energies. We also recall 
that when two partner potentials have continuum spectra the corresponding reflection 
and transmission probabilities are identical. 

The Schrodinger potential ( |4.3| ) that arose in the study of quantum fluctuations on 
domain wall backgrounds has indeed the form V 2 (z) with superpotential 

D - 2 



(i) 



W(z) 



-A'(z) 



(4.15) 



Then, using 
form 



73]) and ( |3.8| ) we find that the partner potential V\{z) in ( [4.1 1|) takes the 



Vi 



f 2/A 
4i? 4 



A 



N 

£ 

i=l 



F-bi 



N 

E 



1(F- 



(4.16) 



which is analogous to the form ( |4.4| ) for V 2 . Specific examples of this relation will be 
considered in detail in later sections. 



5 Algebraic classification 

The underlying mathematical structure for solving the differential equation 

(F'iz))' = (F(z) - b x ){F{z) -b 2 )-.. (F(z) - b N ) , (5.1) 

with arbitrary moduli fej, is that of the Christoffel-Schwarz transformation in complex 
analysis. This transformation is familiar from electrostatics, where one applies the tech- 
nique to find the electric potential for a given distribution of charges. From this point of 
view, it is not surprising that the ansatz we made for constructing static domain walls 
in theories of gauged supergravity amounts to solving a similar mathematical problem. 
It is useful to think of the variable z as being complex, whereas F takes values in the 
complex upper-half plane. Of course, appropriate restrictions have to be made at the end 
in order to ensure the reality of the variable z and hence the reality of our domain wall 
solutions. As we will see in detail, the solutions are characterized by the uniformization 
of Riemann surfaces, which are naturally associated to the Christoffel-Schwarz transfor- 
mation. Hence, the explicit derivation of the Schrodinger potential V(z) requires, for all 
practical purposes, going through such a mathematical framework, apart from its own 
value in providing a systematic classification of all domain-wall solutions in terms of alge- 
braic curves. Also note that the variable F is a function F(zg 2 ), but for simplicity we set 
g — 1 in the following; this parameter can be easily reinstated at the end by appropriate 
scaling in z. Here, we do not assume any particular ordering of 6j. 
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We will treat the Christoffel-Schwarz transformation in a unified way for all three 
cases of interest, namely (D,N) = (4,8) (M2-branes), (D,N) = (5,6) (D3-branes), and 
(D,N) = (7,5) (M5-branes), since there is a hierarchy of algebraic curves within this 
transformation that depends on the isometry groups of the distributions of branes. It 
is useful to start with N = 8 and consider an octagon in the complex z-plane, which is 
mapped onto the upper-half plane via a Christoffel-Schwarz transformation 

dz 

— = (F - bx)-^ lv {F - b 2 )-* 2/n ■ ■ ■ (F — b 8 )'^ /7T . (5.2) 
dF 

This transformation maps the vertices of the octagon to the points bi, b 2 , • • • , b$ on the 
real axis of the complex F-plane, whereas its interior is mapped onto the entire upper- 
half F-plane. The variables tpi denote the exterior (deflection) angles of the octagon 
at the corresponding vertices, which are constrained by geometry to satisfy the relation 
tpi + {p 2 + ■ ■ ■ + (ps = 2tc. We proceed by making the canonical choice of angles <p\ = <p 2 = 
■ ■ ■ = ip 8 = 7r/4, in which case we arrive at the differential equation that relates dz and 
dF: 



(5.3) 



(^) =(F~b 1 )-\F-b 2 )- 1 ...(F-b 8 )- 1 , 

which is the equation we have to solve for the case of M2-branes. 

It is convenient at this point to introduce complex algebraic variables 

x = F(z) , y = F\z) , (5.4) 

which cast the above differential equation into the form of an algebraic curve 

y 4 = {x - bi)(x - 62) • • • (x - b s ) ■ (5.5) 

This defines a Riemann surface of genus g = 9,0 when any two moduli are not equal, 
as follows by direct application of the Riemann-Hurwitz relation, which is standard in 
algebraic geometry. The task now is to uniformize the algebraic curve by finding another 
complex variable, call it u, so that x = x(u) and y = y(u), which resolves the problem 
of multi-valuedness of the algebraic equation above; the corresponding Riemann surface 
is pictured geometrically by gluing four sheets together along their branch cuts. Then, 
following the definition of x and y in terms of F(z) and its z-derivative, one applies the 
chain rule in order to obtain the function z{u) by integration of the resulting first-order 
equation 

dz 1 dx{u) ^ 
du y(u) du 

Finally, by inverting the result we obtain the function u(z), which yields F(z), and 
hence the conformal factor of the corresponding domain wall solutions, as well as the 
Schrodinger potential V(z) for the graviton and scalar field fluctuations in these back- 
grounds. Of course, there is an integration constant that appears in the function z(u), 

5 We will use the symbol g to denote the genus of a Riemann surface. The fact that the same notation 
has already been used for a mass scale should not create a confusion to the reader. 
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but this can be fixed by requiring that the asymptotic behaviour of the domain walls 
approach the AdS geometry as z — > 0. We also note for completeness that there is a 
discrete symmetry x <-> —x, 6j <-> — 6j that leaves invariant the form of the algebraic 
curve. It can be employed in order to set F bigger or equal to the maximum value of 
the moduli hi instead of being smaller or equal to the minimum value, thus insuring that 
z — > corresponds to F — > +00 instead of —00. 

The whole procedure is straightforward, but in practice it turns out to be cumbersome 
when the moduli parameters 6j take general values. After all, the uniformization of a 
genus-9 Riemann surface and the explicit derivation of the function u(z) is a formidable 
task. Matters simplify considerably when one imposes some isometry that effectively 
reduces the genus of the algebraic curve as certain moduli are allowed to coalesce. In 
general we will have models for each continuous subgroup of the maximal isometry group 
SO (8), in which case the algebraic curve takes the irreducible form 

y m = (x - hr (x - hr ■■■(x- b n ) a " , (5.7) 

where the integer exponents (with n < 8) satisfy the relation a\ + a 2 + • • • a n = 2m. 
We present below in table 1 all Riemann surfaces that classify the domain-wall solu- 
tions of four-dimensional gauged supergravity with non-trivial scalar fields in the coset 
SL(8, JR)/ SO(8) by giving their genus according to the Riemann-Hurwitz relation, their 
irreducible form (since in certain cases the exponents have common factors and the curve 
might be reducible when written in its original form), as well as the corresponding isome- 
try groups that determine the geometrical distribution of M2-branes in eleven dimensions. 
We have 22 models in total, which are listed. 
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Genus 


Irreducible Curve 


Isometry Group 


9 


y A = 


(x 


-K 


(x - b 2 ) ■ ■ ■ (x - br)(x - b 8 ) 


None 


7 


y A = 


(x 


-K 


(x-b 2 )---(x- b 6 )(x - b 7 f 


SO(2) 


6 


y 4 = 


(x 


-K 


(x-b 2 )---(x- b 5 )(x - b 6 ) 3 


50(3) 


5 


y' = 


(x 


-K 


■■■(x- b 4 )(x - b 5 ) 2 (x - 6 6 ) 2 


50(2) 


X 


50(2) 




4 


y 4 = 


(x 


-K 


(x - b 2 )(x - b 3 )(x - b 4 ) 2 (x - b 5 ) 3 


50(2) 


X 


50(3) 




3 


y' = 


(x 


-K 


■ ■ ■ (x - b A )(x - 6 5 ) 4 


50(4) 










y 4 = 


(x 


-K 


(x - b 2 )(x - b 3 )(x - b 4 ) 5 


50(5) 










y 4 = 


(x 


-K 


(x — 62)^ — h) 3 (x — b 4 ) 3 


50(3) 


X 


50(3) 






y A = 


(x 


-K 


(x-b 2 )(x-h) 2 (x-b 4 ) 2 (x-b 5 ) 2 


50(2) 


X 


50(2) x 


50(2) 


2 


y A = 


(x 


-K 


(x — b 2 ) 2 (x — b 3 ) 2 (x — 6 4 ) 3 


50(2) 


X 


50(2) x 


50(3) 


1 


y A = 


(x 


-K 


(x - b 2 )(x - 6 3 ) 6 


50(6) 










y 4 = 


(x 


-K 


(x - b 2 )(x - h) 2 (x - 6 4 ) 4 


50(2) 


X 


50(4) 






y 4 = 


(x 


-K 


(x - b 2 ) 2 (x - b 3 ) 5 


50(2) 


X 


50(5) 






y 4 = 


(;r 


-K 


2 (x — b 2 ) 3 (x — 6 3 ) 3 


50(2) 


X 


50(3) x 


50(3) 




y 2 = 


(x 


-K 


(x - b 2 )(x - b 3 )(x - 64) 


50(2) 


1 









y A = 


(x 


-K 


(x - b 2 y 


50(7) 










y = ( 


x - 


-bf 




50(8) 


(maximal) 






y 2 = 


(x 


-K 


(x - b 2 ) 3 


50(2) 


X 


50(6) 






y 4 = 


(x 


-K 


(x - b 2 ) 3 {x - b 3 ) A 


50(3) 


X 


50(4) 






y A = 


(x 


-K 


3 (x — b 2 ) 5 


50(3) 


X 


50(5) 






y = ( 


x - 


-h) 


[x - b 2 ) 


50(4) 


X 


50(4) 






y 2 = 


(x 


-K 


(x - b 2 )(x - b 3 ) 2 


50(2) 


X 


50(2) x 


50(4) 



Table 1: Curves and symmetry groups of domain walls for M2-branes. 



We will see later that, for the models with low genus (0 or 1), the uniformization can 
be carried out in all detail and it is possible in many cases to arrive at explicit expressions 
for the exact form of the Schrodinger potential. 

It is interesting to note that the classification of domain walls of five-dimensional 
gauged supergravity with non-trivial scalar fields in the coset 5L(6, H)/50(6) follows 
immediately from above by restricting our attention to models with an 50(2) isometry 
factor in the symmetry group. It is known that in this case the classification reduces to 
the list of all algebraic curves |l I 



?/ 4 = (x-& 1 )(z-& 2 )..-(z-& 6 ) , (5.8) 

depending on the values of the six moduli 6j. But such curves can be viewed as special 
cases of the N = 8 curves when b 7 = 6 8 = 00. Therefore, by comparison with table 1 we 
list in table 2 all domain walls of five- dimensional gauged supergravity, which correspond 
to various continuous distributions of D3-branes in ten dimensions. We have 11 models 
in total, which are listed. 
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Genus 


Irreducible Curve 


Isometry Group 


7 


y 4 = (x — bi)(x — 


b 2 )---(x- b 5 )(x - b 6 ) 


None 


5 


y 4 = (x — b±)(x — 


b 2 )(x - b 3 )(x - b 4 )(x - b 5 ) 2 


SO (2) 


4 


y 4 = (x — bi)(x — 


b 2 )(x-b 3 )(x-b 4 ) 3 


SO {3) 


3 


y 4 = (x — bi)(x — 


b 2 )(x-b 3 ) 2 (x-b 4 ) 2 


SO '(2) x SO (2) 


2 


y 4 = (x — bi)(x — 


b 2 ) 2 (x - b 3 ) 3 


SO '(2) x 50(3) 


1 


y 4 = (x — b±)(x — 


b 2 )(x-b 3 ) 4 


50(4) 




y 4 = (x — bi ) (x — 


b 2 f 


50(5) 




y 4 = (x — b\) 3 (x - 


-b 2 ) 3 


50(3) x 50(3) 




y 2 = (x — b± ) (x — 


b 2 )(x-b 3 ) 


50(2) x 50(2) x 50(2) 





y 2 = (x — bi)(x — 
y 2 = (x — b) 3 


b 2 ) 2 


50(2) x 50(4) 
50(6) (maximal) 



Table 2: Curves and symmetry groups of domain walls for D3-branes. 



Finally, the algebraic classification of all domain-wall solutions of seven-dimensional 
gauged supergravity with non-trivial scalar fields in the coset 5L(5, H)/50(5) (which by 
the way provides the full scalar sector in this case) follows by considering all Riemann 
surfaces of the form 

y A = ( x -b 1 )(x-b 2 )---(x-b 5 ) , (5.9) 

for various values of the five moduli hi. As before, these surfaces can be viewed as special 
cases of the N = 8 algebraic curves where three of the moduli are taken to infinity, 
i.e. &6 — b-j = bg = oo, whereas the remaining are free to vary. Put differently, we may 
compose the list of all domain walls that correspond to various continuous distributions of 
M5-branes in eleven dimensions by considering all iV = 8 models with a 50(3) isometry 
factor. Thus, we have 7 models in total, which are listed. 



Genus 


Irreducible Curve 


Isometry Group 


6 


y 4 =(x-b 1 )(x- b 2 ) (x - b 3 ) (x - b 4 ) (x - b 5 ) 


None 


4 


V 4 = (x - h)(x - b 2 )(x - b 3 )(x - b A ) 2 


50(2) 


3 


y 4 =(x-b l )(x-b 2 )(x-b 3 ) 3 


50(3) 


2 


y 4 =(x-b l )(x-b 2 ) 2 (x-b 3 ) 2 


50(2) x 50(2) 


1 


y 4 =(x-b l ) 2 (x-b 2 ) 3 


50(2) x 50(3) 





y 4 = (x - b 1 )(x - b 2 ) 4 
y 4 = (x- bf 


50(4) 

50(5) (maximal) 



Table 3: Curves and symmetry groups of domain walls for M5-branes. 



In this latter case, the invariance of the curves under the discrete symmetry x <-> —x, 
h <-> — bi is not present any more because the algebraic equations contain an odd number 
of factors. 



18 



6 Distributions of M2-branes 



In this section we treat the distributions of M2-branes with isometries that correspond 
to genus and 1 and present the uniformization of the associated algebraic curves. We 
have 12 such models of low genus, including the 5*0(8) model of AdS± space. When it 
is possible, we determine the conformal factor of the metrics (written in conformally flat 
form), as well as the Schrodinger potentials corresponding to the equation for the scalar 
field and graviton fluctuations. There are several cases where we can explicitly solve the 
Schrodinger equation and find the spectrum by means of elementary methods, otherwise 
we will use the WKB approximation for the computations. 

We begin first with the model 50 (2) x 50 (2) x SO (2) x SO (2), which is governed 
by a g = 1 algebraic curve, and study certain limits for which some of the moduli are 
let to coincide and the genus is reduced to 0. These cases are SO (4) x 5*0(2) x 50(2), 
50(4) x 50(4) and 50(6) x 50(2) as well as the trivial 50(8) model. We also present 
other g = models, which cannot be obtained as degenerate limits of the above, namely 
the models with symmetry group 50(7) and 50(3) x 50(5). For completeness we also 
present the uniformization of the remaining g = 1 cases, as well as the remaining g = 
model with symmetry 50(3) x 50(4), which unfortunately cannot be explicitly brought 
to a conformally flat frame by expressing their uniformizing parameter as a function of 
z, u(z), in closed form. 



The algebraic curve can be taken from table 1 and corresponds to a g — 1 Riemann 
surface 



with bi = 1/Aj. We set 6 4 = by a shift of x and assume the ordering Ai < A 2 < A 3 < 
A 4 = +oo. Using the birational transformations 



6.1 S0(2) x S0{2) x S0(2) x S0(2) 



y 2 = (x- 1/AiXx - 1/A 2 )(x - 1/X 3 )x , 



(6.1) 



y = Fx 2 /(2v / AiA 2 A 3 ) , X = -l/x , 



(6.2) 



we obtain 



F 2 = 4(X + A!)(X + A 3 )(X + A 3 ) . 
Finally, we can bring the curve to its Weierstrass form 



(6.3) 



w 2 = Av 3 - g 2 v - g 3 , 



(6.4) 



by letting X = v — (Ai + A 2 + A 3 )/3, Y — w. The coefficients turn out to be 



9s = 



92 = 



2 - ((Ax + A 2 - 2A 3 ) 2 + (A 2 + A 3 - 2Ax) 2 + (A 3 + Ax - 2A 2 ) 2 ) , 
4(Ai + A 2 - 2A 3 )(A 2 + A 3 - 2Ai)(A 3 + Ai - 2A 2 ) . 
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This elliptic curve can be uniformized in the standard way using the Weierstrass 
functions^ 

v = p{u) , w = p'{u) . (6.5) 

In general, the two half-periods of the Weierstrass function are given by 

K{k) iK(k') 
ui = , u 2 = , , (6.6) 

V e i - e 3 V e i - e 3 

where K is the complete elliptic integral of the first kind with modulus k and comple- 
mentary modulus k! given as 

k 2 = , k' 2 = 1 - k 2 = e -±^ . (6.7) 

ei - e 3 ex - e 3 

Here ex, e 2 and e 3 are the values of the Weierstrass function at the half-periods, i.e. 
p{uJx) = ex, p(u} 2 ) = tz and p{yj\ + uj 2 ) = e 2 , which are expressed in terms of the 
parameters Aj as 

e i = -A i + ^(Ai + A 2 + A 3 ) , i = 1,2,3 . (6.8) 

So, we finally have 

F(zg 2 ) =x = - ] , (6.9) 

p[u) - (Ai + A 2 + A3J/3 

and 

F'(zg 2 ) = y = ^(AiA 2 A 3 )~ 1/2 p } u) . (6.10) 

2 (p(u) - (Ai + A 2 + A 3 )/3) 2 

Next we compute dx/du = 2(AiA 2 A 3 ) 1 / 2 y and dx/dz = g 2 y and find 

w=y(A 1 A 2 A 3 )- 1 ^ + c, (6.11) 
where the integration constant c is given implicitly by the relation 

p(c) = ^(A! + A 2 + A 3 ) . (6.12) 

The conformal factor is 

^ = —£= 1 1 t [U \ ^^ 2 • (6-13) 
2 V / A7A^ {p{u) - p{c)f 

The constant c in (|6.12j) was chosen in such a way that e 2A ~ 1/z 2 as z — > and the 
space becomes AdS^. Also, according to the conventions followed in this paper we choose 
the branch p'(c) < 0, since c is otherwise determined up to a sign; then, the function F 
in (|6.9| ) indeed approaches +00 as z — > 0. Using the equation for c and certain identities 
involving the Weierstrass function, the conformal factor can be brought to the form 

e 2A = (p(u - c) - p(u + c)) . (6.14) 



6 Throught the rest of the paper we will make use of elliptic functions as well as of other special 
functions following the conventions of EG, O, R2] . 



20 



Finally, the Schrodinger potential ( |4.3| ) can be written as 

V{z) = (2p(« + c) + 2p(u - c) - p(2u)) . (6.15) 

The structure of the spectrum is understood using the results of section 4. For generic 
non-identical values of Aj the spectrum is discrete and the range of < z < z mSuX is finite, 
where z max is found by solving the equation F(z mguX g 2 ) = \/\\. We find, in particular, 
that 

2VA1A2A3 , X , X 

^max = 5 (wi - c) , (6.16) 

9 



where u\ and c are given by (|6.6|) and ( |6.12| ). Then, using ( |4.10| ) with D = 4, A = 4 and 



n = 2, we determine the spectrum within the WKB approximation 

M l = -5—™ ( m + 7j) + £>( m °) , m=l,2,.... (6.17) 

Note that the model we have just considered corresponds to the supersymmetric 
(extremal) limit of the most general rotating M2-brane solution ]24| of eleven- dimensional 



supergravity. Physically, the parameters 6j are set equal (up to a factor) to the four 
different rotational parameters that appear in that general solution, and therefore are 
naturally associated to the Cartan subgroup of 50(8). This identification is analogous 
to the model of a D3-brane distribution with 50(2) x 50(2) x 50(2) isometry, which 
corresponds to the supersymmetric limit of the most general rotating D3-brane solution 



of type-IIB supergravity, as shown in ITT 



We will return again to this model in section 8, together with other potentials that 
arise in elliptic solutions of gauged supergravities, and apply the rules of supersymmetric 
quantum mechanics to simplify the calculation of the exact spectrum. It will turn out that 
the partner potential is formulated as a Lame problem with half-integer characteristic 
n = 1/2. 



6.2 50(4) x S0{2) x S0{2) 

The corresponding algebraic curve of genus is y 2 = (x — bi)(x — 62) (^ — b^) 2 , according to 
table 1. It can be brought to the unicursal form v = w by the birational transformation 

x = hvw ~ 62 , y = (6a - h) v (63 - 63 - (h - b 3 )vw) . (6.18) 
vw — 1 (vw — \y 

Introducing a uniformizing parameter u, so that u = v = w, we find after some compu- 
tation 

2du 

fc = -*,)„■ ■ (6 ' 19) 

which can be easily integrated by elementary functions; the final result depends on the 
sign of (61 - 63) (6 2 - 63). 
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Alternatively, the result can be obtained as the limit Ai = A2 of the SO (2) x 5*0(2) x 
50(2) x 50(2) model. Using some limiting properties of the Weierstrass functions, we 
find that the conformal factor is 



-R^AiA 3 (coslr V^3 - Aiu - A3/A1) 2 

and that the Schrodinger potential takes the trigonometric form 

v - 4^—^- ( ^L Al cosh 2v/- ^ 3 ~ - cosh2 — ~ + 1 - 1 



i? 4 A 2 A 3 V (cosh 2 VA 3 - Xiu - A3/A1) 2 sinh 2 2^A 3 - X lU/ 

(6.21) 

where 

9 2 z sinh'^y/Aa/Ai - 1) 

u = -== + c , c = . (6.22) 

2A1VA3 v A3 — Ai 

The spectrum is continuous (according to the ordering of our parameters Ai = A2 < A3), 
and < z < 00. This example corresponds to n = 4, which we have already discussed 
on general grounds in section 4. Hence, there is a finite mass gap given by the value of 
the potential at z = 00. We find, in particular, that 

ML = 4fe^ • (6.23) 



A 2 A 3 i? 4 



Alternatively, one may consider the limit A2 = A3 of the model with 50(2) x 50(2) x 
50(2) x 50(2) symmetry. It turns out that the result is equivalently described by 
considering the analytic continuation of ( |6.21| ) and ( |6.25| ) by taking A3 < Ai. Then, the 
conformal factor becomes 



e 



2.4 



(A 3 /A 1 ) 1/2 (l - A3/AO 3 / 2 , (6.24) 



R 2 Xi\ 3 K a/ LJ y " LJ (cos 2 VAi - A 3 n - A3/AO 2 
and the Schrodinger potential is 



i? 4 A 2 A 3 V (cos 2 VAi - X 3 u - A3/A1) 2 sin 2 2VAi - A 3 w, 

(6.25) 

where 

g 2 z shr^y/i - A3/A1) 

u = = + c , c = . (6.26) 

2AiV A3 v Ai — A3 

The spectrum is discrete and the range of z is < z < ;z max , where 



2A lv / A7sin- 1 ( v /A 3 7Ai") 



(6.27) 



gWXi - A 3 

Using the general results found in (|4.10|) (with D = 4, A = 4, n = 2), we obtain within 
the WKB approximation the spectrum 

M l = m ( m + tt! + 0( m °) , m=l,2,.... (6.28) 



2 ■■- y - ■ 2 

max x _ 
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The exact solution : It might seem hard to find the explicit solution of the Schrodinger 
equation with potentials ( |6.21| ) and ( |6.25| ). However, we will show that the explicit 
solution can be readily obtained by applying techniques of supersymmetric quantum 
mechanics, thus going beyond the WKB approximation. We consider for simplicity only 
the case when the spectrum is discrete and the potential is given by ( |6.25| ). It turns out 
that the partner potential to ( |6.25| ) is given by 

„ 4(Ai - A3) ( 3 \ 

using the notation of section 4, where u is given by ( p.2t ). If we change variables for the 
wave-functions of the potential V\ to 

x = cos 2 a/Ai — A3U , < x < ^ < 1 , 

gr(i) = x 3/4 (1 _ X f/^Y{ X ) , (6.30) 

the new independent variable Y (x) satisfies the hypergeometric equation with parameters 
a = 3/2 + yU, 6 = 3/2 — and c = 2 in the standard notation, where the constant /i is 
defined as 

Hence, we may write down the solution for the wave-function that is regular at 

x = 0, as 

~ x^(l - xf'F Q + ^- - K 2; x\ . (6.32) 

Demanding that it vanishes at the end point located at x = A3/A1 < 1, yields the 
quantization condition for the spectrum 

3 A \ 

+ Mm, 2 ~ Hm, 2; y J = , m=l,2, .... (6.33) 

It is not possible to find a closed formula for the exact spectrum except for the particular 
case when Ai = 2A3. Then, we have < x < | and making use of the identity 



3 3 1\ _ Vtt 

2+^2" /im ' 2; 2 J " r(| + ^)r(f-^) ' (6 ' 34) 

we see that the condition ( |6.33| ) (with A3/A1 = 1/2) is satisfied provided that [i m = 2m+|, 
where m = 1, 2, . . .. Finally, by employing ( |6.31| ), we find that the mass spectrum is given 
by 

16 1 

A 1 = 2A 3 : M 2 m = ^m{m + -) ) m = l,2,.... (6.35) 

One may easily check that it coincides with the WKB result (|6.28|) (after setting 
Ai = 2A3), which is also exact in this case. The situation is similar to the distribution of 
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D3-branes with SO (2) x 50(4) symmetry for which the exact spectrum || |9| coincides 
with what has been found using the WKB approximation ||25|| . This is a characteristic 
property of shape-invariant potentials j|8] . We also note presently that the solution 
( f3.32| ) can be expressed in terms of Jacobi polynomials as 



Ai =2A 



3 • 



< x < - 
- - 2 



m 



1,2,... . (6.36) 



For the generic case, where the constant parameters Ai and A3 satisfy the condition 
Ai > A3, but otherwise they remain unrelated, one may solve ( |6.33|) in the asymptotic 
regime of large \i. Using the fact that for large \i 



3 A 3 



r(/z - 1/2) 



y u- 1/2 (sin^)- 3/2 sin/i^(l + C( / u- 1 )) , (6.37) 



-3/2 



r(/i + i/2)' 

where the angle ip is determined by cos ip = 1 — 2A3/A1, we see that the spectrum is given 
to leading order by the solutions of the simple trigonometric equation sin fup = 0. It 
follows easily that the result agrees with the leading term in the WKB formulae ( |6.28|) . 

(2) 

Finally, note that the wave-functions \l/m of the corresponding partner potential 
V2 = V have a form that is more complicated than ( |6.32| ). They are obtained using 
( |1.13| ), where for the lowering operator ( 4.12 ) we have W(z) = A'(z) with A(z) given by 
(fT24j). The result is 



X 



X 



Aire - A 3 



8 (A 3 - 2A 3 x + Aix 2 ) F 



2 "I" /V> 2 ^ tm ' 



2,x + 



+ (Afi 2 m - 9) (x - x 2 ) {X x x - A 3 ) F ( b - + /i m , b - 



(6.38) 



with x being defined in ( |6.30|) ; at the two endpoints W m vanishes 



6.3 50(4) x 50(4) 

This model is described by an algebraic curve with g = and can easily be obtained by 
taking the limit A3 — > 00 in the above expressions. We find that the conformal factor for 
the metric is 

, (6.39) 



whereas the potential is 



In this case we have < z < 00, and the spectrum is continuous with a mass gap 

4 

i? T A 2 
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Alternatively, the same result is obtained by direct uniformization of the irreducible 
genus curve, y — (x — bx)(x — b 2 ) (see table 1). Letting v = y/(x — bi), w = x — b 2 , we 
arrive at the standard form v = w, which is uniformized by a complex parameter u as 



x 



u 



b 2 , y = u(u + b 2 -b 1 ) 



and so the final result reads 
z(u) = - 



In 



u 



u + b 2 - bi 



u(z) 



bi 



(6.42) 



(6.43) 



choosing the integration constant to be zero. Then, we have 



y 



(h - b 2 f 



4sinh 2 |(6i - b 2 )z 



(6.44) 



which is equivalent to the expression for the conformal factor above, after introducing 
the appropriate scale. 

In this model we can also determine the eigenf unctions exactly. It turns out that the 
wave-functions are expressed in terms of hypergeometric functions, using a characteristic 
parameter q ^ 0, as 



where 



y 



C 2 y l - q F 



cosh 



1 q q . I- 



ClV 1+q F I -1-1 -1,1- q] 



+ 



1 q q 1 l 
"2 + 2'2' 1 + g; ? 



2z 



Q 



M 2 R 4 Xl 



(6.45) 



(6.46) 



' ' V ' i 

with constant coefficients C\ and C 2 . A solution valid for q = can also be written down 



following |5l| , but it will not be needed for the present purposes. Because of the mass 
gap (|6.41 ), the parameter q is purely imaginary. This provides an orthonormalizability 
in the Dirac sense (with the use of a 5-function), since the solution ( 6.45| ) becomes an 
incoming and an outgoing plane wave for z — > oo. On the other hand, examining the 
behaviour of ( |6.45| ) near z = 0, we require the coefficient of the most singular term to 
vanish. It yields the following condition, 



Ci 



r(l-g) 



Co 



r(i + ?) 



r (I - 1) r (i - 1) -r(| + |)r(i + f) 



o 



(6.47) 



Note that the supersymmetric partner of the potential ( |6.40|) is just a constant (given 



by V, 



and therefore, according to the results of supersymmetric quantum me- 



chanics, the potential ( |6.40| ) is reflectionless jiS 
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6.4 S0(2) x 50(6) 



In this case the irreducible algebraic curve is given by 

y 2 = (x - bi)(x - b 2 ) 3 . (6.48) 

Then, the birational transformation 

wvb 2 -b 1 Oi-M 2 (RA a\ 
x = 7- , y = -w- , 6.49 

wv — l {wv — \y 

brings it to the standard unicursal form v — w, which is uniformized as usual by in- 
troducing a complex parameter u, u = v = w. We may easily find the coordinate that 
brings the four-dimensional metric to a conformally flat form, 



2(u - 1) 
Z -- g z{b x -b 2 ) 

with conformal factor given by 



(6.50) 



(2 - g 2 (h - b 2 ) z 



e 2A = v g v 1 L > ' 2 . (6.51) 
g 2 z 2 (4 — g 2 (bi — b 2 ) z) 

In turn, the Schrodinger potential becomes 

V(z) = l 9 4 (h~b 2 ) 2 2gHh-b 2 ) 2 

z 2 4 (2 -g 2 (6i-6 2 ) ^) 2 (4-^(^i-& 2 )^) 2 ' 

For &! < 6 2 we have < z < oo and the spectrum is continuous with no gap, which is in 
agreement with our general discussion. Otherwise, for b 2 < b±, we find < z < z max = 
2/ (g 2 (bi — b 2 )) and the spectrum is discrete. Using (|4.10 ) with D = 4, A = 4 and n — 2, 



we find within the WKB approximation that the spectrum is 



Ml = h)2 m f m +^\+ Q(m°) , m = l,2,... . (6.53) 



Note for completeness that Q6.51 ) and ( |6.52 ) can also be obtained from (|6.20 ) and ( |6.21 ) 



and from (|6.24|) and (|6.25| ) by simply taking the limit A3 — > Ai. 

Finally, let us mention that the wave-functions and the spectrum can be determined 
exactly using, as before, techniques of supersymmetric quantum mechanics. The same 
is true for the other two cases that will be discussed in subsections 6.5 and 6.6 below. 
Since there are certain similarities in all three cases of interest, we choose to present 
their analysis based on supersymmetric quantum mechanics all together, and uniformly, 
in subsection 6.7. 
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6.5 50(3) x 50(5) 

The algebraic curve of this model with genus is given by 

y 4 = ( x - hfix - b 2 f , (6.54) 
which by means of the following birational transformation 

w 3 vb 1 -b 2 4 {b 1 -b 2 f ( . 

x = — = — , y = -w v— , (6.55) 

w a v — i (w 6 v — \y 

can be brought to the standard form w = v , which is uniformized (as usual) by a 
parameter u, u = v = w. The coordinate choice that brings the four-dimensional metric 
to conformally flat form is 

4(u - 1) 



g 2 (&i -b 2 )u ' 
and the corresponding conformal factor is 



(6.56) 



e» = 1024(4-^ 

^ 2 z 2 (8 -^ 2 ^) 2 (32- g 2 bz (8-g 2 bz)f 



with b = bi — b 2 . 

The Schrodinger potential turns out to be 

1 ^! ( 3 8 
^ ~~ 4 V(4-(? 2 6z) 2 + (8-(7 2 6^) 2 " 

512 16 



(32-g 2 bz (8-g 2 bz)Y 32-g 2 bz (8-g 2 bz) 



(6.58) 



For b\ < b 2 the spectrum is continuous without a gap, in which case we have < z < oo. 
For bi > b 2 , we find < z < z max = A/(g 2 b) and the spectrum is discrete; within the 
WKB approximation, it is given by the simple expression 

M ™ = 1^™ ( m + I) + ° (m0) ' rn^l,2,... , (6.59) 
where (|4.10| ) has been used with D = 4, A = 4, n = 3. 



6.6 50(7) 

The corresponding g = algebraic curve is given by 

y 4 = ( x - b^ix - b 2 ) (6.60) 



and by means of the birational transformation 

w 3 vbi-b 2 (h-b 2 ) 2 

x = — = — , y = -w— , 6.61 

w 6 v — 1 [w 6 v — l) z 
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it can be brought to the unicursal form v = w , which is again uniformized as u = v = w. 
The coordinate that brings the four-dimensional metric to conformally fiat form is 



4(1 -u 3 ) 

V(6i-&2) 



(6.62) 



and the conformal factor is given by 



2l fpfi-W*)* 



e 2A = — * . (6.63) 

(-l + (l-f^^)tj 

The potential of the Schrodinger equation of this model is 

g A b 2 (-5 + 42 (1 - \g 2 bz) (l - \g 2 bz)* +91 (l - \g 2 bz)~A 

V(z) = ^ — '- , (6.64) 

64 (-1 + \g 2 bz+ (1 -\g 2 bz)^ 

with b = b 2 — b\. For b\ > b 2 we have < z < oo and the spectrum is continuous without 
a mass gap. On the other hand, for b\ < b 2 , we have < z < z m3iX = 4/(3g 2 b) and the 
spectrum is discrete given, within the WKB approximation, by 

Ml = ^^m(m + ^] +O(m ) , m = l,2,... , (6.65) 

where ( |4.10| ) has been used with D = 4, A = 4, n = 1. 



6.7 Some exact results 

We are in a position to apply techniques of supersymmetric quantum mechanics, as in 
subsection 6.2 above, in order to find the explicit solution of the Schrodinger equation 
with potentials ( |6.52| ), ( |6.58| ) and ( |6.64| ) and determine the associated spectra in a uniform 
way; all three quantum mechanical problems will be treated at once. 

First consider the case where the spectrum is discrete. Let us change the variable to 
x, z — z max -i, where z max is defined as in the appropriate subsections and < x < z max . 
Then, by the rules of supersymmetric quantum mechanics, the partner potential to our 
problem is given by 

z/ 2 -l/4 

Viix) = , (6.66) 

x z 

where v = 1,2, ±2/3 for the three models that correspond to SO (2) x 5*0(6), £0(3) x 
50(5) and 50(7), respectively. Both signs ±2/3 were chosen for the 50(7) model for 
reasons that will be indicated towards the end of this subsection. The solution of the 
Schrodinger equation which is regular at x = (for v = 1,2,2/3) is given in terms of 
Bessel functions as 

~ x 1/2 J v (Mx) . (6.67) 
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For v = —2/3 the wave-function diverges at x — 0, but it is not strong enough to make 
it non-integrable.[] Imposing the condition that the wave-function vanishes at x — z max , 
which also ensures the Hermiticity of the Hamiltonian, we find the mass spectrum in 
terms of the zeros of the Bessel function, 

Ju(M m z max ) = , m = 1,2,... . (6.68) 



Then, the states in ( |6.67|) , which correspond to the different solutions of (|6.68| ), constitute 



a complete set of states. The undetermined overall constant in ( |6.67|) can be found, as 



usual, by demanding the orthonormalizability condition dx^n^m = & n ,m and 

using the fact that 

dxxJ v {M m x)J v (M n x) = ^J,+i(M m z max )<5 m , n . (6.69) 



Note that the conditions (|6.68| ) and v > — 1 are crucial for the validity of this equation. 



An asymptotic expression for the eigenvalues M m can be found, which is valid for 
large values of the argument M m z mssL of J u . Using standard formulae from the theory of 
Bessel functions we arrive at the result 

M 2 m = JL— m (m + v - + O(m ) , m=l,2, .... (6.70) 

2 max V 2 / 



For v — 1 and v — 2, it agrees with the WKB formulae (|6.53| ) and (|6.59 ). However, the 



WKB formula (|6.65|) is reproduced for the value v = —2/3 (by first shifting m by one 
unit) instead of v — 2/3. Actually, for v = 2/3, we obtain (|6.65| ) with the number 5/6 
replaced by 1/6. This ambiguity in the spectrum remains unresolved even for the partner 
wave-functions that we are interested in computing afterall; these are obtained using 
( |4.13j ), where for the lowering operator in (|4 . 1 2|) we have W(z) = A'(z), with A(z) given 



by ( |6.51| ), ( |6.57|) and (|6.63| ) for all three different cases respectively. 
We find, in particular, for the wave-functions that 

(2x 
2 _ 2 Ji{M m x) + M m J {M m x) 
Z max X 



V 



( 4r 3 \ 
~ 4 _ M M m x) + M m J x {M m x) , (6.71) 

V ^max J 



v = ±\- ~ ±if 2 _ 4/3 (x 4/3 J± 5 /3(M m x) + z^J Tl/3 (M m x)) 

2max X 



It is crucial to note here that for all three type of wave-functions above, vanish 
at the end point x = z max - This can be shown by first expanding the wave-function 



7 In fact, the general criteria developed in ]53| | render the propagation of quantum test particles in 
such a space-time geometry as unphysical, and hence should be related to unphysical vacuum expectation 
values of the scalar fields. Indeed, it was shown in |0] that, in this case, the density of the M2-brane 
distribution has a negative component which is physically unacceptable. 
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around x = z max , using properties of the Bessel functions, and then observe that the 
coefficient of the divergent part in the expansion is proportional to Jv(M m z mgcK ), and 



hence vanishes due to (|6.68|) ; the constant part vanishes identically. Note also that the 
asymptotic behaviour of the wave-function near x = is ~ x 1 ' 6 for both v = 2/3 
and v = —2/3. Hence, there is a priori no reason to dismiss either one of the two values 
z/ = 2/3orz/ = —2/3. It remains unclear to us, at least for the moment being, what is the 
extra condition one should impose in order to exclude one of these two values. As soon 
as this becomes possible, the spectrum corresponding to the distribution of M2-branes 
with SO (7) symmetry will be determined unambiguously in the discrete case. 

Considering the same potentials ( |6.52| ), ( |6.58[ ) and fl6.64p , but for continuous spec- 
trum, we find in all three cases that the parameter z max becomes negative. By changing 
variable to x, where z = x — |-2 m ax|, we note that the appropriate partner potential is 
still given by (|6.66|) , but with \z max \ < x < oo. Since the point x = is not contained in 



the range of x we should admit both independent solutions of the corresponding Bessel 
equation, namely J v and N v , following the standard nonmenclature. Demanding that the 



wave-functions vanish at the endpoint x 
Hence, the wave-function is given by 



determines their relative coefficient. 



~ N u {M\z max \)J u {Mx) - J V {M\ z max \)N„(Mx) 
whereas the spectrum is continuous with no mass gap. 



(6.72) 



6.8 50(3) x 50(4) 

Last, but not least, we consider the remaining curve of genus 

y 4 = (x-& 1 )(x-6 2 ) 3 (z-& 3 ) 4 , (6.73) 
which can be brought to the unicursal form using the birational transformation 

b\ — b 2 v 3 w {b\ — b 2 )v fb\ — b 2 v 3 w 



x 



V 3 W 



V 3 W 



1 — V 3 W 



(6.74) 



and hence uniformized by setting u = v = w. Assuming for definiteness that b\ ^ b 2 , we 
may proceed to solve the differential equation (|5.6|) in order to obtain the corresponding 
function z{u). Taking the limits b\ — > b 3 or b 2 — > b 3 yields the models £0(3) x SO (5) 
and SO (7) respectively, which we have already discussed. For general frj's we arrive at 



a < 



1 



1 



2 tan" 1 



V2q 



u 



+ ln 



\[2qu + q 2 u 2 



a > 



1 



9 2 (bi - b 3 )q 



1 + \f2qu + q 2 u 2 

—2 tan -1 qu + In 



q 2 u 2 



+ const. 



+ 



(6.75) 



(r 


+ qu\ 




- qu) 



+ const. 
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with 



a = — , q = I ct 1 4 . (6.76) 

h - h 



Unfortunately, it is not possible to invert the relations and find u(z) in closed form, 
and so they will not be pursued any further. 

6.9 More g = 1 surfaces 

In the following we present the uniformization of the remaining models with genus 1 in 
a unifying way. These cases are 

(i) y 4 = (x- a )(x-b){x-cf : 50(6), 

[ii) y A = {x-a){x-b){x-c) 2 {x-d) A : S0{2) x 50(4) , 

(Hi) y A = (x-a)(x-cf(x-bf : 50(2) x 50(5) , 

(iv) y A =(x-c) 2 (x-af(x-bf : 50(2) x 50(3) x 50(3) . 

Using birational transformations, they can be brought into the same form 

(X - cfY A = (X - a) (X - b) , (6.77) 

where in each case we consider the following: 

(i) X = x , Y !J 

(ii) X = x , Y 
(Hi) X = x , Y 



(iv) X = x , Y 



(x 


-c) 2 ' 






y 




(x 


— c)(x 


-d) 




y 




(x 


— b)(x 


-c) 


(x 


— a)(x 


~b) 


y 



Then, the birational transformations are employed 

v 2 | 2a6-c(a+ 6 ) £ 
~y 2 c(a—b) 16 

= ° 7 .2 I k a+b-2c , fc2 ' 
U 2 a-b U T 16 



b — cw , 9 (a — 6) 2 (a — c) 

y = fes- w,th fe = - (t-cy • (6 ' 78) 



to bring the common form (|6.77 ) into the Weierstrass normal form of the curve 



k 2 

w 2 = Av 3 - g 2 v- g 3 ; g 2 = — , ^3 = . (6.79) 



31 



Note that in all cases we have set x — X, with X given by (|6.78 ) above, whereas for 
y we have to treat each transformation separately. Explicit calculation shows that y is 
equal to 

" 3 — 1 



(0 



(ii) 



[III] 



[IV] 



[a — c 



vw 



,2 ' 



y/(a-c)(b-c) 2 + |o±^ + »y 

(a - c? w( a ~ d > + (c~d) J^ c (v + k/Af 



v / (a _ c)(& _ c)2 ^2 + ka±^c v + ^y 

same as in (ii) setting d = b , 

(a - c) 2 2v(u + A;/4) 2 ( a ~ 6 ) u + ( c ~ b )\f^( v + W 



(b 



w 



( 2 i k a+b-2c , k±\ 



As usual, in all four models we have to solve the differential equation (|5.6j ) in order to 
determine the function u(z). 

Using v = p(u) and w = p'(u) we find the following results in each case separately: 

4 



a — c 



CM + ^ 



\(b -c) ; 



+ const. 



(c - d)(a - 6) 2 ^/(a-^^-c) . 



-k 2 u+ 



V + — V- 

v+p'(a-) 

V + — V- 

4 



<j{u — a+) . , 

log _ _. , _ + { +2{(a + )u\- 



log 



er(u + a + ) 
oiu — d-\ 



(6.80) 



<j\u + a_ 



+ 2C(a_)u 



+ const. 



a — b 
2 



2((u) + 



PW 



p(u) - k/A 



+ const. 



-u + const. , 



where v± are the two roots of the equation 



(v + k/A) 2 + 



a — d a — c 







(6.81) 



c — d\ b — c 

and a± are defined by v± = p(a±). In the case (ii) above we assume that c ^ d, 
and generically a± differ from the half-periods of the corresponding Riemann surface. 
Although case (i) corresponds to taking c = d, there is no smooth limit of (ii) that yields 
the expression (i). In the case (iii) above we have v+ — and so the derivation has to 
be performed separately without taking (ii) in the limit b = d. 

It is rather unfortunate that we can not invert the relations and find u(z) in closed 
form for these models, apart from the case (iv). Hence, they will not be discussed any 
more. We only leave case (iv) as an exercise for the interested reader to explore it further. 
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7 Distributions of M5-branes 



In this section we treat the distributions of M5-branes with isometries that correspond 
to genus 1 and and present the uniformization of the associated algebraic curves in 
as much the same way as for the M2-branes. There are only two models to consider 
apart from the AdS 7 space, namely the distribution with isometry group 50(2) x 50(3) 
that corresponds to a genus 1 algebraic curve, and that with isometry group 50(4) 
corresponding to a genus algebraic curve. Unfortunately there is no case where we 
could find the exact spectrum in terms of known functions and therefore we will only 
resort to the WKB approximation for the computations. 



7.1 50(2) x 50(3) 

Consider the genus 1 curve of table 3 

y 4 = (x-a) 2 (x-b) 3 . (7.1) 

With the aid of the birational transformation 

b — aw 2 . w ( w 2 \ , 

x = ^-r +a ' y = ^ b - a K{^- 1 ) ' (7 - 2) 

it can be brought into the Weierstrass form 

w 2 = Av 3 - , (7.3) 

which is uniformized using the Weierstrass functions defined for g<± = and g% = 0. 
Since 

x = b + ( Q ~ 6 ) 2 = F{z ft (a- b) 2 p'(u) = , 2 

X ° + 16p(w) 2 " [ 9 J ' V 16p(u)3 " [ 9 j ' 1 ] 

we obtain 

m = c-^, (7.5) 

where c is an integration constant (to be fixed by the asymptotic conditions). Then, 
using (|3.3| ), we find the conformal factor 



> 2A = g 2 F\zg 2 )l . (7.6) 



Note that the two half-periods of the torus U\, u 2 can be computed using ( |6.6j ) and 
the fact that in our model the modulus and its complement are equal to each other, 
k = k' = l/y/2, since g 3 = 0. We have 

r(i/4) 2 

Ux = - iuJ2 = ±Stm . (7.7) 



Hence, for a > b, u)\ is real and 102 is purely imaginary, whereas for a < b they are complex 
conjugate of each other. The constant of integration c is determined by requiring that 
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the space becomes AdS 7 for z — » 0. We find, in particular, that c — u± + u 2 , which is 
complex (or real) if a > b (or a < b). 

We also find the following Schrodinger potential, using the rescaling factor g — 2/R, 



(7-8) 

Note that < z < z max , where z max = ^- + u 2 ) if a < b and z max = ^-oj\ if a > b. 
The spectrum is discrete and, within the WKB approximation, it is given by 

167r 3 y/b — a 
r(l/4) 4 i? 4 



a<b - Mm^-—i——m(m + 3) + O{m ), m = l,2,..., (7.9) 



a>b: Ml = ^J^rn(m + 2) + 0(m°) , m = l,2,... , (7.10) 

where ( [4.1U| ) has been used with D = 7, A = 4 and n = 3 (or n = 2) for a < b (or 
a > b). We will see later, in the context of supersymmetric quantum mechanics, that the 
partner potential is related to the potential of the 5*0(3) x 5*0(3) model of D3-branes in 
five dimensions. 



7.2 50(4) 

Consider next the genus curve 



y 4 = (x — a)(x — b) 4 . 



The birational transformation 

1 



x 



V 3 W 



+ a 



y = - 



+ a — b 



(7.11) 



(7.12) 



V \v°w 

brings it into the unicursal form v — w, which can be uniformized with a complex 
parameter u, as v = w = u. Consequently, we arrive at 

4 du 



dz 



g 2 (a - b)u 4 + 1 
which yields upon integration the following cases: 



(7.13) 



a = b 
a <b 

a > b 



4it 

-— + const. , which gives AdS 7 , 
T 



In 



g 2 yb — a \_ \ 1 — u\Jb — a 



+ 2tan" 1 (wv / ^ r 



+ const. ,(7.14) 



g 2 V2^a~^b 
+ 2 tan" 1 



In 



1 + uV2^a~^b + u 2 <fa^b 
1 - uVztfa^b + u 2 <fa^b 

uV2^a~^b 



+ 



1 — u 2 \/a — b 



+ const. . 
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Apart from the maximally symmetric model that corresponds to a = b, these relations 
cannot be inverted to yield u{z) in closed form. 

The nature of the spectrum depends crucially on the sign of a — b. Using our general 
formulae we find that for a < b we have < z < oo and that the spectrum is continuous 
with a mass gap given by 



a < b : M| ap = . (7.15) 



For a > b we have < z < z max = \/2ir/g 2 (a — 6) 1 / 4 and the spectrum is discrete. It is 
approximated by the WKB formulae ( |4. 10| ) with D = 7, A = 4 and n = 1 as 



a>b: M 2 m = 8v ^7^ m(m + \) + O(m ) , m = l,2,.... (7.16) 
Finally we would like to mention the relation of the models we have presented in 



this section to the most general solution of rotating M5-brane |26| of eleven-dimensional 
supergravity. The latter, besides the usual Poincare invariance along the brane, has also 
an 50 (2) x 50 (2) symmetry group corresponding to the Cartan subgroup of 5*0(5). 
Hence, in the extremal limit it will correspond to a supersymmetric solution associated 
with an algebraic curve of genus 2, as can be seen from the appropriate entry in table 
3. The two independent parameters in the equation of the algebraic curve are related 
to the rotational parameters of this rotating M5-brane solution. The genus 1 model 
with symmetry S0(2) x 50(3) corresponds to the particular limit when one rotational 
parameter is set equal to zero. The case with a < b corresponds to the rotating solution 
with Lorentzian signature, whereas for a > b it corresponds to the same solution, but 
with the time and angular parameters analytically continued so that the metric remains 
real but its signature becomes Euclidean. The associated spectra are given in the WKB 
approximation by ( |779|) and ( pMO|) fl The genus model with symmetry group 5*0(4) 



corresponds to letting the two angular parameters become equal. We note that this is not 
equivalent to setting one of them zero and keeping the other finite. Then, a < b describes 
the Euclidean solution, whereas a > b describes the Lorentzian. The corresponding 
spectra are described by ([7.1 5|) and (|7.16|) respectively. 



7.3 Wilson surfaces 

We would like to calculate the vacuum expectation values of Wilson surface operators 
in the six-dimensional (0, 2) theories on the Coulomb branch. It was shown in that 



the AdS/CFT correspondence could be used to compute Wilson surface observables |jB| 
of (0, 2) theories in the limit of a large number Nj, of M5-branes. The Wilson area 
operator in the supergravity picture is defined by requiring that a membrane ends at 



8 In fact ( 7.10 ) coincides with the supersymmetric limit of the WKB formula given in |5J] for the 
masses of ++ glueballs using rotating M5-branes with one rotational parameter, in a supergravity 
approach to QCD4 B . 
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the boundary of AdS-j x 5 4 on the surface that defines the operator. We will consider 
Wilson surfaces corresponding to a pair of parallel strings on the boundary using the 
prescription of [55] in the special backgrounds constructed in sections 3.1, 7.1 and 7.2. 
Wilson loops turned out to be useful tools for learning about the physics of gauge theories 
in the study of supergravity duals of four- dimensional theories on the Coulomb branch. 
It is interesting that complete screening was found with an associated screening length 
suppressed by l/^g s Nb compared to what is expected from field theory considerations 
at weak coupling [g, [|. 

In the conformal limit, this calculation was performed in |55|] leading to the result 

- ~F(ir t? • (7 - 17) 

where E denotes the energy per unit length or tension between the infinitely long strings 
as function of their separation L. This can easily be generalized to backgrounds of the 
form ( |3.15| ), ( |3.18| )- fl3~20| ) with / defined as in equation (|3.5|). We will impose the minor 
restriction h\ = . . . = b n = b and b n+ i = . . . = 65 = 0, thus breaking the 5*0(5) symmetry 
to SO(n) x 5*0(5 — n). Furthermore, we choose the orientation of the Wilson surface on 
the deformed 5 4 to be constant and to lie in the subspace spanned by y n +i, ■ ■ ■ , 2/5 (see 
also (I3TT8D ). 

Minimizing the membrane action in these backgrounds with the orientation chosen 
as above, yields the following integrals for the length and the energy (for 1 < n < 4) 

L = * , ^™ , (7.18) 

h = (F-6)?fV , g = F' 1 , (7.20) 

with F Q > b being the minimal value of F that the Wilson surface reaches. In general these 
integrals cannot be expressed in terms of known functions and so we will only present here 
some numerical results. For F ^> b, or equivalently for small separations L R 2 / \/b, 
the behaviour of the potential is as in (|7.17f ) and goes to zero faster for larger separations. 
From a certain distance L > L max , and further on, there does not exist a minimal surface 
connecting the two strings on the boundary. Instead, a configuration of two separated 
surfaces hanging straight into the interior of the geometry is energetically preferred. This 
means that the potential is screened for large separation; a phenomenon that was also 
observed in four-dimensional super conformal theories on the Coulomb branch |8], The 
maximal distance at which the string breaks can be determined numerically 

R 2 

Lmax C n (7.21) 

with ci ~ 0.71 , C2 ~ 0.69 , C3 ~ 0.70 and C4 ~ 0.78. For n — 1, 2, 3 the length reaches 
its maximum L max at F > b and becomes zero as F — » b. At L = L max the energy 
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is larger than zero and the split configuration is preferred. For n — 4, L max is reached 
exactly when F — * 6; at this point the potential tends smoothly to zero and will remain 
there even if the separation is increased. This is the phenomenon of complete screening 
that was also found in the context of some special continuous distributions of D3-branes 
in |8], |9|. This phenomenon occurs in cases where the mass spectrum is continuous with 
a mass gap. 



8 Comments on Lame equations 

In this section we summarize some results on the Lame equation, and its various gener- 
alizations, which arise in the study of quantum fluctuations for the scalar and graviton 
fields in the background of domain walls associated to elliptic functions. We have already 
seen that in many cases the Schrodinger potential has the common form 

V(u) = \{X + l)p{u)+^(fI+l)p(u + UJ 1 )+u(u + l)p(u + U 1 + U2)+K{K+l)p{u + UJ 2 ) , (8.1) 

for various choices of the coefficients A, fi, v and k; they are all constrained to satisfy the 
Hermiticity bound > —1/2. We have found, in particular, the following list of examples: 

(z) D3 -branes with SO (2) x SO(2) x SO{2) : A = |, fi = v = k = -- , 

3 1 

(ii) D3 — branes with 5*0(3) x £0(3) : A = v — - , /i = K = - , 

5 1 1 

(izi) M5-branes with S0(2) x S0(3) : A= 2' V= 2~' ^= K =~^- 



The first two cases were derived in [fTI |, where emphasis was placed on analyzing domain 
wall solutions of five- dimensional gauged supergravity.^] There, the expressions for the 
conformal factor were found to be 



6 -{2W) Wp(u)J ' (8 ' 2) 

with u = z/R 2 and u = z/(2R 2 ), respectively. The third example was discussed here in 
section 7.1, and has u = 2z/R 2 . The parameter u assumes real values from to uj\ (real 
semi-period) for the cases (i), (ii) and (iii) (with a > b), whereas for the case (iii) (with 
a < b) u takes real values from to u\ + u>2- In the context of supersymmetric quantum 



mechanics j|8] the potential (i) is mapped to a potential with coefficients A = 1/2, 
fi = v = k = 0, which is simpler to study. The potentials (ii) and (iii) (with a > b) also 
turn out to be related to each other via supersymmetric quantum mechanics, and the 
details are worth exposing. 



9 Actually, for the SO(3) x 50(3) model of D3-branes, the Schrodinger potential was originally pre- 
sented in another form in fll|| , but that is equivalent to the potential (ii) because of special identities of 
the underlying Riemann surface with 173 = 0. 
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More precisely, we find that the supersymmetric partner potentials corresponding to 
the SO (2) x 5*0(2) x SO (2) model in D = 5 are given respectively by the pair 

Vl(u) = 3p(2 M ) , V 2 (u) = Ap(u) - p{2u) , (8.3) 

making use of the identity 

= ^(p(u) + ptu + u^ + p(u + u 2 ) + p(u + u 1 + u 2 )) . (8.4) 

So, by rescaling -u by a factor of 2, setting -u = 2u, we obtain a partner Schrodinger 
problem in u with potential V\{u) = A(A + l)p(u) having A = 1/2 as advertised. For 
later use we drop the tilde and still use the variable u, but with a range from to 2u>i. 
According to the previous general discussion, the ground state energy of Hi is expected 
to be zero, which will also be encountered later using a direct approach. 

The Schrodinger potentials appearing in the examples (ii) and (iii) (with a > b) also 
have supersymmetric partners within the same class. We find, in particular, that the 
supersymmetric partner potential of the model (ii) has A = v = 1/2, /i = k = 3/2, 
whereas a similar analysis for the model (iii) yields a potential with A = v = 3/2, 
H — k — 1/2. We do not observe important simplifications occuring in the form of the 
potential in any of the two cases. An interesting observation is that the supersymmetric 
partner of the potential (iii) is (ii), and we believe that there is a deeper reason for this 
finding. Note at this point that, generically, any two supersymmetric partner potentials 
are related to each other by simply changing sign in the superpotential W <-> — W . For 
domain wall solutions this would mean that the conformal factor reverses, since A <-> —A, 
and so this transformation could only be of mathematical interest in relating different 
spectra. In physical terms, the transformation W <-> —W cannot be used to map one 
domain wall solution to another because it fails to preserve the AdS boundary condition 
imposed on the conformal factor for z — > 0. However, it may happen in certain cases (as 
above) that there are two different superpotentials with the correct asymptotic behaviour 
as z — > which yield the same supersymmetric partners; indeed, the potential (ii) is the 
V\ partner of a potential V 2 given by the model (iii), due to special identities on Riemann 
surfaces with g 3 = 0, but conversely this is not so because W <-> — W does not relate the 
model (ii) and (iii). 

Another notable relation concerns the supersymmetric partner of the elliptic potential 
for the SO (2) x SO (2) x SO (2) x SO (2) model of M2-branes in four-dimensional gauged 
supergravity. Actually, in this case we find that the two partner potentials 

Vi = 3p(2u) , V 2 = 2p(u + c) + 2p(u - c) - p(2u) (8.5) 

are connected by supersymmetry, and so by reinstating the overall scaling factor and 
constant shift that relates the uniformizing parameter to the Schrodinger variable (now 
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called u), we obtain again the Lame potential with A = 1/2 but in the range from 2c 
to 2o>i. This concludes the presentation of some qualitative results on the Schrodinger 
equation of quantum fluctuations on elliptic backgrounds. 

An interesting problem that remains unsolved is the exact evaluation of the full spec- 
trum of the Schrodinger equation in this class of potentials. So far we have relied on 
semi-classical approximation methods to get a feeling about the spectrum and the exis- 
tence of a mass gap. Unfortunately, the exact result is very hard to find, even in some 
simple cases, and involves transcendental equations in rather implicit form. Nevertheless, 
it is quite instructive to highlight some special results in order to appreciate the degree 
of difficulty one faces in the general case. We note that the family of potentials under 
investigation are indeed a natural generalization of Lame's potential n(n + l)p{u) by 
adding terms located at all four corners of the parallelogram u, u+U\, u + u>2, u + u)\ + ui2 
in the complex domain of a genus 1 Riemann surface. The Lame potential was originally 
introduced in order to describe the analogue of spherical harmonics for the solutions of 
the Laplace equation in three dimensions with ellipsoidal symmetry. As such, n can only 
take (positive) integer values. However, generalizations were also considered, more than 
a century ago, for half-integer and other values of the parameter n. The main results 
in this direction go back to Hermite, Brioschi and Halphen (see, for instance, [|57j and 
f58fl), but also Darboux apparently studied some aspects of the general potential ( |8.1|) in 
its Jacobi form. Half-integer values of the coefficients are particularly interesting for the 
examples we have at hand, but we will be able to say something explicit only for the case 
A = 1/2 and y. = v = k = 0, which is related to the SO{2) x SO(2) x SO{2) model of 
five-dimensional gauged supergravity or to the 5*0(2) x 50(2) x 50(2) x 50(2) model 
in four dimensions. 

We proceed by considering the Lame equation 

d 2 



th{1 + n(n + l)p(u)J = EV(u) , (8.6) 

where E are the energy levels of the corresponding one-dimensional quantum mechanical 
problem and ^(u) are normalized wave functions that vanish (typically) at u = modulo 
the real period 2u>i. Of course, when n is a positive integer, the solutions are easily 
described using appropriate ratios of the Weierstrass sigma-function, namely 

,a(u - Pi)cr(u - (3 2 ) ■■■a(u- (3 n ) 



au 



^i{u) = e 

(T"[U) 

^2{u) = e — , (8.7) 

<j n [u) 

where the constants a and $ are determined by substituting the ansatz into Lame's 
equation. Using these two (in general independent) solutions one can construct regular 
solutions by taking suitable linear combinations of them and obtain a transcendental 
equation for the energy eigenvalues. 

On the other hand, if n is not an integer the solutions will be difficult to describe, 
even in a formal sense. Some simplifications occur when n is half of an odd positive 
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integer. It is known in this case, using the substitution 

= (p'(z))~ n $(z) , where z = u/2 , (8.8) 
that Lame's equation transforms into the differential equation 

§- 2 "^|f + 4K2„-l)p W + B )* = 0. (8.9) 

Then, according to results obtained by Brioschi and Halphen more than a century ago, 
a formal solution can be written as 

00 

^z) = J2c r (p(z)-e 2 ) a - r , (8.10) 

r=0 

provided that 

(a-2n)(a-n + l/2) =0 , (8.11) 
and that the following recursive relations are satisfied: 

(a-r-2n)(a-r-n+l/2)c r + (3e 2 (a - r + l)(a - r - 2n + 1) + e 2 n(2n - 1) + c r _! 

= (ei - e 2 )(e 2 - e 3 )(a - r + 2)(a - r - n + 3/2)c r _ 2 . (8.12) 

Here, ei, e 2 and e 3 denote the three roots of the cubic curve in its Weierstrass form. 
Actually, when n is half of an odd positive integer, there is a solution expressible in finite 
form 

71- 1/2 

<&(*) = £ c r (p(*) - e 2 ) 2 "^ , (8.13) 

r=0 

which corresponds to a = 2n, and provides discrete energy levels E by solving the 
recursive relations with c n+ \/ 2 = 0; they are all real for curves having real roots e^. 
Otherwise, for solutions expressible as an infinite sum, the energy levels remain arbitrary 
by these general considerations alone. 

The simplest case to consider has n — 1/2. It turns out that a solution expressible in 
finite form has only one term: 

$(z) = Co(p(z)-e 2 ) , (8.14) 

with energy E = and Co arbitrary. This yields one of the two independent solutions 
of Lame's equation with zero energy, call it ^i(u), whereas the other is obtained by 
employing the general formula that relates any two solutions with the same energy E: 



where C is another constant. An explicit calculation shows that the most general solution 
with E = is given in terms of two integration constants A, B by 

*(«) = -4= (Ap(u/2) + B) . (8.16) 
VPW 2 ) 
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Having established this, one may impose the regularity of the physical solution at u — 0, 
namely ^(O) = 0, and set the coefficient A = 0; otherwise the wave function will diverge 
as l/y/u for u — > 0. For u = 2cui, however, the wave-function blows up and hence it is 
not normalized. 

For n = 1/2, but with E ^ 0, the formal solution has an infinite number of terms. 
According to the general discussion, we find after some calculation (working with either 
a = or a = 1) the result 

= C (± + y - , (8.17) 

where C is an arbitrary constant and 
di = -g, ^ = -(^- + e 2 

ds = ~ ((ei-e 2 )(e 2 -e 3 )+9(f + e 2 ] ( ^ + e 2 ) ) , (8.U 



12 v V 6 7 V 18 

d t = l(( ei -e > )(e > - e ,)g + 2e,) + 9(| + e ,)g + e 1 )(| + ej 

and so on. Unfortunately, even in this simplest case with n = 1/2, it is very difficult to 
solve the recursive relations and find all coefficients dk in closed form in order to sum 
up the infinite series of terms. In any case, this procedure yields one formal solution of 
Lame's equation with E ^ 0, ^i(u), whereas the other (independent) formal solution 
of the same energy, ^(w), can be obtained according to the general formula above. 
However, the integration that determines \£ / 2 (m) in terms of ^i(u) cannot be performed 
explicitly unless one knows first how to sum up the infinite series of terms for &(z). 

Note that for n = 1/2 the formal solution ^i(u) appears to be regular at u — 0, 
whereas at u — 2uj 1 (where p{u/2) = e±) this is not guaranteed. In fact, since the 
coefficients dk are polynomial functions of the energy E, one must demand that the 
resulting "energy series" converge. Since p'(oui) = 0, we demand i.e. 

(8 - i9) 

This is certainly a non-trivial constraint on the allowed energy bands when E ^ 0, which 
also depend on the relative size of the a- and 6-cycles of the Riemann surface, i.e. the 
differences e 2 — e 3 and e 2 — e\ that appear in the "energy sum" . However, we are not able 
at present to find the complete solution to the problem in closed (even transcendental) 
form. 



9 Conclusions 

We have investigated in detail the structure of domain wall solutions in theories of gauged 
supergravity in diverse dimensions by considering the effect of non-trivial scalar fields tak- 
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ing values in the coset space SL(N, JR)/ SO(N). The presentation was kept quite general 
to cover, where possible, aspects of domain wall solutions in theories of D- dimensional 
gravity by turning off the effect of any other fields, such as gauge fields and fermions. 
Special emphasis has been placed on two cases, namely (D, N) = (4, 8) and (7, 5), which 
arise by compactification of eleven-dimensional supergravity on S 7 and S A respectively. 
The effect of the scalar fields in four (or seven) dimensions is related to deformations of 
the round spheres in the compactifying space, thus breaking the isometry group SO (8) 
(or 5*0(5)) into appropriate symmetry subgroups. In fact, we were able to give an alge- 
braic classification of all such cases using the Christoffel-Schwarz transformation, which 
arises in the solution of the first-order Bogomol'nyi-type equations for the conformal fac- 
tor and the scalar fields of these models. As a result, for D = 4, we found a hierarchy 
of 22 solutions starting from an algebraic curve of genus 9 corresponding to completely 
broken isometry group. When some cycles shrink to zero size, by letting some moduli 
coalesce, the symmetry group is enhanced, whereas the genus of the Riemann surface is 
lowered accordingly. Among the genus models there is the AdS 4 space with no scalar 
fields, having maximal isometry group SO (8), while in all other cases AdS 4 is reached 
only asymptotically. Similarly, for D = 7, we found 7 models in total, which can be 
classified starting from a genus 6 algebraic curve with no isometry and proceeding all the 
way down to genus 0, where the AdS 7 space arises with maximal isometry group SO (5). 
All other models in the list admit subgroups of SO (5) as isometries and approach AdS-j 
only asymptotically due to the presence of non-trivial scalar fields. 

A geometrical picture of our solutions in four and seven dimensions is provided in 
terms of distributions of M2- and M5-branes in eleven dimensions. The analysis has 
been carried out in detail choosing suitable harmonic functions that describe continuous 
distributions of branes in eleven dimensions, and is in accordance with the geometric 
deformation of seven and four dimensional spheres induced by the non-trivial moduli 
of the underlying algebraic curves. The resulting picture resembles the construction of 
domain walls of five- dimensional gauged supergravity with non-trivial scalar fields in the 
coset SL(6,1R)/SO(6)i where suitable continuous distributions of D3-branes were con- 
sidered in ten-dimensional type-IIB supergravity. An interesting problem that remains 
open for further study is the possible effect of dualities on the structure of domain-wall 
(and other) solutions in various dimensions. For example, M2- and M5-branes are related 
to each other via electric-magnetic duality, and hence various distributions of them in 
eleven dimensions should yield strong-weak coupling relations among solutions of gauged 
super gravities in lower dimensions. Also, D3-branes appear in a sequence of dualities be- 
tween extended objects in higher dimensions, using both S- and T-dualities, and so the 
variety of domain walls in D = 4, 5 and 7 dimensions (where consistent truncations of 
supergravity are known to exist) ought to be interelated. Of course, it will be interesting 
to formulate this in the algebro-geometric context provided by the Christoffel-Schwarz 
transformation, where there is a universal moduli space of the genus 9 Riemann surface 
for M2-branes y A = (x — bi)(x — 62) • • • (x — b$). It has been noted earlier that the Riemann 
surfaces of D3-branes arise as special cases by taking two of the moduli to infinity, while 
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for M5-branes there are three moduli taken to infinity. Thus, it is quite natural to expect 
that the modular transformations of the underlying algebraic curves have a natural inter- 
pretation in higher dimensions as S- and T-dualities among branes and their continuous 
distributions thereof. We hope to return to this issue in a separate publication. 

Another topic that has been investigated in detail concerns the analytic form of fluc- 
tuations for the graviton and scalar fields in our domain-wall backgrounds. Since the 
quantitative analysis of this problem amounts to solving a one-dimensional quantum me- 
chanical problem with potential V(z) that behaves as 1/z 2 for z — > (the asymptotic 
AdS region), the right identification of the variable z and the associated Schrodinger 
potential V(z) become crucial for extracting the spectrum. We found that the variable z 
appears naturally in the Christoffel-Schwarz transformation, which is a complex trans- 
formation that maps the interior of a closed polygon in the z-plane onto the upper-half 
F-plane, while the variable F is more appropriate for the brane description of our con- 
figurations in higher dimensions. Thus, the uniformization of the associated algebraic 
curves is a necessary step in order to derive the exact form of the potential, in each case 
of interest, and, consequently, to determine its spectrum. We were able to complete this 
mathematical task for the models with Riemann surfaces of genus and 1 and derive 
the potentials in closed form for most of these cases. Subsequently, using techniques of 
supersymmetric quantum mechanics, which ensure that the spectrum is non-negative, we 
were able to compute the spectrum exactly for many models and make estimates using 
the WKB approximation for many others. The elliptic models exhibit an interesting 
class of potentials within the family of generalized Lame potentials, which can be studied 
analytically only in certain cases. 

It might seem surprising that we have encountered domain wall solutions with very 
little or even no isometry in the classification in terms of Riemann surfaces. Of course, we 
are unable to perform the uniformization and find explicit expressions for the conformal 
factor of their metrics and their scalar fields in general, since this has only been done 
for models of low genus and hence bigger isometry groups. Nevertheless, this situation 
can be perfectly accommodated in supersymmetric theories, where all the domain-wall 
solutions leave half of the supersymmetries unbroken. It is rather instructive to compare 
this situation with the geometry of four-dimensional hyper-Kahler manifolds and their 
possible isometry groups. Recall that all hyper-Kahler manifolds are supersymmetric 
(preserving half of the supersymmetries) and as such they admit three independent com- 
plex structures I, J, K, hence a whole sphere of them, since al + bJ + cK will be a 
complex structure if a 2 + b 2 + c 2 = 1. The group 5*0(3) acts naturally on the space of 
complex structures as rotations; but this does not necessarily imply that all hyper-Kahler 
manifolds are 5*0(3) symmetric. Although many examples (like Eguchi-Hanson, Taub- 
NUT and Atiyah-Hitchin) have an 50(3) isometry group, there are others with less 
isometry or none; for example Dancer's manifold has only an 50(2) isometry, whereas 
K3 commonly used in compactifications of string theory has no isometries at all. In 
this sense, our domain-wall solutions to gauged supergravities in four, five and seven 
dimensions provide useful tools for developing a deeper understanding of the consistent 
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truncations of eleven- dimensional supergravity, and they can be rather exotic. Also, the 
associated Bogomol'nyi bounds and their possible description using contour integrals on 
the underlying Riemann surfaces pose some interesting mathematical questions for the 
future. 

Finally, there is the conceptual issue of relating 7?.-symmetry to the isometry group 
that remains unbroken by the geometric structure of our solutions, which we belief is 
worth emphasizing. In our approach, the ^-symmetry is spontaneously broken by giving 
vacuum expectation values to the scalar fields of the theory that are charged under this 
symmetry. Then, the 7£-symmetry is not a symmetry anymore, but relates different 
vacua of the theory. Generically, this procedure also breaks conformal invariance, thus 
breaking half of the 32 supersymmetries independently of the amount of symmetry that is 
left unbroken. In supergravity this can be understood as follows: for generic values of the 
real parameters 6, in the master curve y A = (x—bi)(x—b2) ■ ■ ■ (x—bs), the gauge symmetry 
(related to the 7£-symmetry in field theory), is spontaneously broken. In fact, for any 
given choice of the parameters, the curve is not invariant under SO (8). However, its form 
is preserved since SO (8) acts naturally on b%, b%, ■ ■ ■ , bs in its fundamental representation 
and rotates any given choice of moduli bi into another. This is different from the situtation 
where the background describes a flow from the maximally supersymmetric theory to a 
conformal theory with less supersymmetry. In the latter case the theory is perturbed 
by adding suitable deformations to the action, which explicitly break some or all of 
supersymmetries; at the IR fixed point the geometry is again AdSo and the 7£-symmetry 
is related to the number of supersymmetries. This means, in particular, that the isometry 
group contains a factor that is equal to the 7£-symmetry of the field theory. 
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